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This  thesis  describes  a  scene  analysis  algorithm  which 
locates  targets  in  a  noisy  background.  Preprocessing  is 
used  to  detect  the  edges  of  objects  in  a  digitized  scene. 
The  edges  extracted  fros  the  scene  is  then  cross-correlated 
with  a  teaplate  of  the  target  to  be  found. 

Cross-correlation  ie  done  by  using  , cesplex-ccn Jugate 
siultlplication  in  the  frequency  doaaln.  Areas  of  high 
correlation  are  recerrelated  using  snaller  isuges  which  can 
ee  processed  faster.  The  potential  targets  are  energy 
norsalised  with  respect  to  the  teaplate  in  order  to 
elialnate  false  correlation  with  noise. 

Bven  in  a  scene  with  high  energy  noise,  the  algorltha 
works  very  well  when  operated  under  the  constraints  of  size, 
orientation,  and  perspective.  /'  .'h>\  )r\ T<Ti '  -.~U <?<r! 
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SCENE  ANALYSIS 


USING  RECUESIVE  FREQUENCY  DOMAIN  CORRELATION 
WITH  ENERGY  NORMALIZATION 

I .  INTRODUCTION 

1*1  Background 

Iaage  pattern  recognition  ie  concerned  with  extracting 
diicrialnaMng  features  froa  visual  data  and  using  these 
features  to  identify  objects  in  the  iaagc.  The  huaan  aind 
can  perfora  this  function  with  little  trouble;  however, 
aachines  cannot.  The  applications  of  iaage  peitern  recogni¬ 
tion  lie  priaarily  in  four  aajor  areas:  (Ref  •) 

1 .  Docuaent  processing 

2.  Industrial  autoaation 

3.  Medicine  and  biology 

A.  Military  target  finding/guidance* 

This  thesis  deals  aalnly  with  the  ailltary  applications 
of  target  recognition;  however,  the  results  presented  in 
tnie  thesis  can  be  extended  to  the  other  ereas. 

With  recent  advances  in  visual  and  infrared  detector 
technology,  the  problea  of  collecting  iaage  data  Is  a  rela¬ 
tively  ainor  obstacle.  The  aain  problea  is  the  interpreta¬ 
tion  of  an  iaage:  filtering  out  the  unwanted  visual  data  to 
deteralne  if  a  target  1  r.  present. 

Although  there  has  been  auch  work  done  in  this  area,  no 
perfect  algoritha  has  yet  been  produced  which  can  process 
the  visual  data  collected.  One  reason  is  became  the  level 


of  the  scene  clutter  (noise  energy)  varies  for  every  scene. 
This  thesis  deals  siainly  with  the  elimination  of  all  noise. 

Much  work  has  bean  done  on  this  area  at  the  AFIT  Signal 


Processing  Laboratory.  This  thesis  is  a  continuation  of 
four  recent  theses  done  at  AFIT:  Horev,  Haaadonl ,  Croaer, 

and  Wells  (*ef  3,2,1,10);  as  well  as  the  proposed  disserta¬ 
tion  by  McKeeaan  (Ref  7). 

1.2  ftoblts 

An  iaaga  of  a  scene  that  contains  a  potential  target 
has  aore  inf oraation  that  just  the:  fron  the  target.  For 
exaaple,  infrared  radiation  is  detected  fros  the  targets,  as 
well  as  froa  the  background  of  the  scene.  The  coaputer 
which  handles  the  visual  data  aust  be  able  to  distinguish 
the  target  froa  the  background  rad*  .tlon. 

The  problea  of  locating  targets  in  a  scene  has  six 
aajor  variables: 

1.  Target  position 

2.  Target  shape 

3.  Target  slse 

A.  Target  perspective 

3.  Target  rotation  (orientation) 

6.  Scene  clutter  (noise). 

The  deteralnat ion  of  all  target  positions  within  a  scene  is 
the  goal  of  this  thesis.  The  target  shape  is  the  property 
which  characterises  the  target.  The  siie  and  perspective  of 
the  target  are  also  variables  within  the  problea,  but  they 


■ay  be  known  In  a  |lvan  acanarlo  (for  example,  obttrvlnt 
tanka  froa  an  airplane  of  known  altitude).  The  aain  thruat 
of  thia  atudy  waa  to  locate  the  targeta  regardleee  of  the 
aaount  of  ecene  clutter  (nolee). 

The  atandard  aethod  for  the  collection  of  thia  viaual 
data  la  to  repreaent  tha  acena  with  an  array  of  pixala 
(picture  eleaenta).  Tha  etandurd  iaage  file  used  at  tha  APIT 
Signal  Proceaaing  I.ab  containe  65,536  pixala  (256x256). 
Each  pixel  containe  A  bite.  Thia  allowe  for  16  different 
gray  lavala  (0  to  15)  which  corraapond  to  tha  intanaity  of 
the  radiation  at  each  pixel  point. 

When  the  iaaga  ia  taken  at  a  relatively  cloaa  dlatanca 
froa  the  acane  tha  targeta  aay  be  dietinguiehable  to  a  huaan 
obeerver.  When  the  iaage  ia  taken  at  a  relatively  far 
diatan'-e  froa  tha  acena  tha  targeta  appear  mm  blobe.  Often 
tlaea  even  a  huaan  obeerver  cannot  identify  tha  blobe  aa 
being  targeta. 

Tha  problaa  ia  to  develop  an  algorithm  which  aanlpu- 
latee  the  iaage  data  eo  that  the  targeta  aay  he  located  end 
identified  with  little  or  no  huaan  intervention* 

1.3  Scope 

Thia  theela  deal#  only  with  the  proceaaing  of  the  iaage 
data  collected  froa  a  acena.  The  aethod  of  collection, 
optica  Involved,  ataoapherle  tranaaiaalon ,  and  wavelengtha 
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The  scop*  of  this  thtiit  confines  tho  eoftwere  to  be 


used  on  the  computers  within  the  AFXT  Signal  Processing  Lab. 
The  elgorlthns  were  tested  on  vleuel  image  dote  created 
within  the  lab,  although  the  process  la  not  neccessar i ly, 


confined  to  visual  data. 


in  a  two- 


dimensional  array  of  pixels,  visual  or  Infrared,  can  be 
operated  on  by  these  algorithms. 

Constraints  on  the  computer  run  times  wars  not  consi¬ 


dered  in  this  thssis< 


iasge  processing  slgorlthm  which  achieved  the  goal  of  target 
identification.  The  philosophy  taken  was  that  if  the  proper 
algorithm  could  be  developed  in  the  Signal  Processing  Lab, 
regardless  of  run  time,  then  it  may  ha  Implemented  in  a  fast 
running,  dedicated  VLSI  system  (haf  4). 

The  software  written  has  the  flexibility  to  wort  with 
several  sixes  ef  video  data  files.  However,  due  to  computer 
run  times  with  the  larger  files,  many  of  the  secondary  steps 
use  only  files  sixes  of  44x44  pixels. 


1.4  Assumptions 

The  algorithms  developed  perform  under  the  following 
aseumptlons  concerning  the  input  test  scenes i 


1.  T» 


iltlpli 


Igitii 


gets  which  were  similar  in  slse,  persective,  and  orienta¬ 


tion. 


2.  The  slse,  perspective  and  orientation  of  the  tar* 
gets  contained  within  the  scans  is  known.  As  stated 
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earlier,  size  and  perspective  may  be  given  in  a  specific 
scenario.  The  problem  of  orientation  may  be  solved  by  re¬ 
analyzing  the  rcene  for  each  possible  orientation  of  the 
target.  This  is  believed  to  be  much  the  same  way  that  the 
human  mind  prc  . eases  images  that  it  sees  (Ref  4). 

3.  The  location  and  relative  energy  of  the  targets  is 
not  known.  The  amount  of  scene  clutter  is  also  unknown. 

4.  When  only  images  of  size  64x64  were'  used  for  recor¬ 
relation  it  is  assumed  that  the  target  was  no  more  that  64 
pixels  in  size  in  any  direction  or  that  an  area  on  the 
target  of  that  size  contained  the  neccessary  identifying 
visual  information. 

5.  The  unknown  locations  of  the  targets  may  be  any¬ 
where  in  the  scene.  However,  Indent  1 f icat ion  of  targets 
close  to  the  edge  of  the  scene  may  be  somewhat  difficult 
compared  to  those  which  are  at  least  -Keif  the  distance  of 
the  recorrelation  file  size  away  form  the  border. 

6.  The  digitized  scene  images  may  contain  noise  (scene 
clutter).  This  noise  ,nay  be  of  high  or  low  energy  level. 

The  separation  of  individual  targets  is  great 
enough  to  allow  for  only  one  target  to  be  contained  in  each 
of  the  smaller  images  used  for  recorrelation. 


v 
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The  proctii  of  target  recognition  described  in  thie 
theeie  is  divided  into  four  aajor  parts: 

1.  Preprocessing  of  the  scene  to  eliminate  background 
noiee  and  enhance  the  target  edges.  The  edges  of  the  ob¬ 
jects  are  used  to  constrain  the  identifying  characteristics 
of  the  potential  targets  and  to  reduce  overall  scene  clut¬ 
ter.  Edge  detection  id  also  used  becauee  it  is  a  forn  of 
low-level  processing  perforaed  by  the  hunan  visual  eystea. 

2.  First  pass  correlation  of  the  edge-detected  scene 
with  a  teaplate  of  the  target  to  be  located.  This  correla¬ 
tion  is  done  to  identify  potential  target  locations  and  le 
perforaed  in  the  frequency  doaain. 

3.  Second  pass  correlation  (in  the  frequency  doaain)  of 
potential  target  locations  with  the  target  teaplate.  Energy 
noraallsetion  is  perforaed  on  the  potential  target  to 
enhance  correlation  with  true  targets  and  attenuate  correla¬ 
tion  with  false  targets. 

4.  Evaluetior.  of  the  recorrelat ion  video  files, 
l.ft  Egu i ynent 

All  equipaent  used  in  this  theeie  ree»5ee  in  the  API? 
Signal  Precaeeing  Laboratory.  All  software  written  seaains 
on  aagnetic  tape  at  the  lab.  Appendis  A  contains  a  brief 
description  of  the  equipaent  and  software  support  aveilable 
in  the  lab. 
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II 


PREPOCESSIMG 


2.1  Video  File* 

Most  of  the  loftmre  in  this  thesis  is  designed  to 
operate  on  video  files  of  the  sis*  ISitlSi  pixels*  Each 
pixel  contains  an  integer  value  corresponding  to  the  Inten¬ 
sity  of  the  radiation  at  that  particular  point  In  the  scene. 
These  values  rang*  fros  0,  for  total  black*  to  13*  for  total 
whit*.  The  software  nay  be  aodified  to  handle  ioages  of 
different  sixes  or  of  different  pixel  resolution*. 

This  pixel  resolution  requires  4  bite  to  represent  each 


point.  For  econosy  of  storage  the  video  files  were 
condensed  so  that  each  14-bit  word  contains  the  infornatlon 


of  4  pixels.  The  video  files  were  configured  in  this  aanner 
with  the  subroutine  REPACK  and  retrieved  with  the  aubroutine 
UNPACK.  The  source  cod*  for  these  subroutine*  is  Contained 
in  Appendix  D. 

2.2  Edge  Detection  Algorlthu 

before  the  targets  sre  extracted  freu  the  scene*  a 
preprocessing  algorlthu  is  perforsed  to  reduce  noise.  The 
Wedge  function  wo*  used  to  do  this  by  detecting  the  visible 


was  found  to  be  the  nost  effective  and  the  sost  consistent 
of  the  edge  detection  operators  that  have  been  studied  (Kef 
7).  The  Wedge  algorlthu  used  was  that  of  Wells's  work  in 
(tef  to). 


\ 


7 


.  .  p  V.  ^  *»  S  -*►  W  , 


This  process  has  two  stops.  Tho  first  Is  the  extrac¬ 
tion  o  f  edges  In  f  our  directions.  The  second  step  is  the 
combining  of  the  four  edge  extractions  into  one  video  file. 

The  first  step  is  performed  by  the  program  WEDGE.  This 
program  accepts  the  name  of  a  video  file  at  d  the  names  of  4 
output  files  to  be  created.  Each  of  tho  four  output  fllae 
contains  the  edges  of  the  scene  in  a  certain  direction 
(vertical,  horliontal ,  and  tho  two  diagonale).  WEDGE  uses 
several  3x3  maeks  to  detect  edges  of  pixel  intensity.  This 
is  implemented  by  the  subroutine  WE0W1,  which  usee  a  four 
element  wedge  operator  (Eef  10).  figure  2.1  shows  tho  masks 
used  by  WE0W1 .  Each  of  tho  four  directions  has  four  masks 
which  are  ueed  to  detect  edges  in  the  given  direction. 

The  program  WEDGE  examines  t!.e  value  of  each  plxal  and 
the  pixels  which  surround  it.  Each  of  tha  nine  pixel  values 
are  multiplied  by  the  corresponding  coefficient  of  the  mask. 
Into  one  of  thp  four  output  files  is  put  a  value  which 


indicates  tha  level  of  tho  greatest  degree  of  sharpness  of 
tho  edge  in  tho  corresponding  direction.  This,  procoes  is 
continued  for  every  pixel  point. 


Figure  2.2  contains  a  234x234  pixel  representation  of  a 
group  of  F-14'e.  Figures  2. 3-2. 4  contain  the  four  output 


video  fllae  of  the  four  different  masking  operations. 


•>V: 

V-„” 


Direction  #1  (H< 


»i) 


@ 


4~- 

-♦---4 

4-— 4  — 

—  4—  —  —  4 

f  2 

J  2 

i  5  i 
i  *  i 

I  -  *  !  -l 

}-l  1 

-4-- 

-4---4 

4---4«- 

-4-— —4 

1-1 

1-1 

1-1  1 

l-i  i-i 

!-»  1 

a— 

-4— -4 

4---4-« 

—4— ——4 

i-1 

l-l 

1-1  1 

t  2  |  2 

t  2  | 

-4---4 

4—4-- 

—4— ——4 

Direction  #2  i 

(let  Diagonal) 

-4-— 4 

4— -4--. 

-4—4 

I-I 

1  2 

J  2  | 

1-1  1-1 

|-»  1 

♦ — 

4-— 4--. 

-4— ——4 

1-1 

1-1 

J  2  J 

1  2  J-l 

|-»  1 

4-~. 

-4—4 

4-—  4— 

*4—4 

1-1 

1-1 

1-1  1 

t  2  t  2 

1-1  1 

♦ - 

»4 - 

-4 - 4 

-4 - 4 

Direction  #3  (Vertical) 

•  4 - 4 

4—4— 

J-l 

1-1 

1  2  J 

J  2  J-l 

1-1  1 

4 - 

•4-— 4 

4—4— 

•4—4 

1-1 

1-1 

1  2  J 

1  2  J-l 

1-1  1 

4--- 

•♦—  ——4 

•4—4 

1-1 

1-1 

1  2  | 

1  2  |-1 

J-l  1 

4--- 

•4"—  — 

•4— —  —4 

4 4— 

4-— 4 

♦  - 4 - 4 - 4 

{-2  5-2  |-2  { 

♦  - 4- — 4—4 

»  »  J  1  I  1  J 

♦- — ♦ - ♦ - 4 

J  »  i  I  I  l  J 

♦  - 4- —4 — -4 


—  —  -4- 

.--4- 

- 4 

1  J 

l  ! 

1  I 

- 4- 

•  --4 

1  ( 

1  1 

1  1 

•  —4* 

•  --4 

1-2  |-2  |-2  | 


t  I  I  1  I  I  | 
4— 4— 4— 4 
1-2  |  I  |  1  | 

4— 4— -4 - + 

i-2  |-2  |  1  | 

— 4---4 


»  »  I  l  1-2  J 

4— 4— 4— 4 
I  »  I  1  J-2  | 

♦ - 4 - 4 - ♦ 

J  1  J  1  J-2  J 

4— 4«— «4— 4 


1-2  |  I  |  I  | 


£ 


i 

\ 

* 

3 

j 

‘3 


i 

I 


r. 

r. 

r. 

r 


Direction  #4  (2nd  Diagonal) 


4  — ;*■ 

-4— 

-4... 

-4 

.4--- 

-4 - 4 - 4 

4- 

... 

t-1 

J-l 

J-l 

1 

1  2 

1  2  |-1  J 

l 

1 

4--- 

-4 - 

-4 

4--- 

♦ 

i 

i 

i 

4 

1 

1 

1 

♦ 

a 

4- 

»  •  * 

1-1 

J-l 

J  2 

1 

1  2 

J-l  J-l  l 

1 

1 

4—  —  > 

*4— — « 

*♦- - « 

-4 

• 

♦ 

i 

i 

i 

4 

l 

l 

l 

♦ 

4- 

... 

J-l 

t  2 

i  2 

J 

J-l 

J-l  J-l  J 

J 

1 

4--- 

■4—  — « 

•4 

4--. 

-4—4—4 

4- 

1-2  J-2  J  1  | 

l“2  i  1  f  1  J 
4— -4— 4— 4 

I  »  I  »  I  I  I 


Figure  2.1 
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Figure  2.2  256  x  256  Pixel  Representation  of  F-I6*s. 
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The  iccita,  containing  P-lia,  in  Figure  2.2  ia  tho 
•xaaplt  uicd  throughout  tho  taut  of  this  thaaiv.  This  iaage 
vaa  used  for  aevtral  reasons : 

1.  It  contained  Multiple  targets  which  net  the  pre¬ 
viously  stated  constraints  in  the  Assuaptiona. 

2.  It  had  a  noisy  background. 

3.  It  contained  high  energy  noise  froaa  the  words  arti¬ 
ficially  placed  on  the  scene  (the  scene  was  originally  a 
■agazine  cover). 

The  second  step  of  the  edging  algoritha  is  the  coa- 
bining  of  the  four  output  files  into  one  output  file  which 
contains  the  edges  in  all  directions.  This  is  aceoapllshed 
by  the  prograa  COM*  which  eoabines  the  four  output  files  of 
the  progran  WEDGE.  Figure  2.7  show  the  result  of  COMB 
perforwed  on  Figures  2. 3-2. 6. 

The  software  for  the  Wedge  preprocessing  is  ,  contained 


in  Appendix  B 
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*  £i**l  **§§  CO|||LATION  PROCESS 
3.1  TEMPLATE  GENERATION 

(•fort  Che  correlation  proccit  be|in«  a  tceplate  of  the 
deaired  target  auit  first  be  generated.  The  teeplatei  were 
first  taken  directly  froe  the  acene  iteelf  and  later  froa 
other  aourcea.  In  any  caae,  the  teaplate  will  exiat  in  aoae 
video  file.  Proa  this  video  file  aoae  optional  preprocea- 
aing  using  the  Octek  board  aay  be  used.  The  teaplate  auat 
be  in  a  rectangular  area  that  is  free  of  any  unwanted  noise. 
The  prograa  OCTEK  has  the  ability  to  block  out  and  replace 
any  area  of  the  video  file  with  any  other  gray  level  (in 
this  case  0  for  black).  Pigure  3.1  above  the  deaired  tea¬ 
plate  inf oraation  froa  Pigure  2.7  after  the  noiae  was 
reaoved  using  OCTEK. 

The  teaplate  video  file  .is  produced  by  the  prograa 
SPIN.  Thia  prograa  creates  a  new  video  file  which  places 
the  teaplate  at  the  origin  of  the  coordinate  eyetea.  In  the 
frequency  doaain  the  teaplate  aeswaeo  to  exist  in  a  wrap¬ 
around-space.  Thia  naans  if  an  iaage  aovea  off  the  screen 
at  the  top,  it  will  reappear  at  the  bottoa;  if  it  aovea  off 
the  screen  to  the  left,  it  will  reappear  on  the  right,  at 
cetera* 

Using  this  concept,  SPIN  places  the  original  teaplate 
at  the  origin  of  the  new  acene.  Thia  will  result  in  the 
correlation  peaks  having  the  saae  coordinates  of  their  cor- 
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responding  targets  or  noise.  This  is  done  by  quartering  the 
original  tcaplate  and  aoving  the  four  portions  to  tha 
opposite  corners  in  the  neu  teaplate  scene.  The  area  con¬ 
taining  the  teaplate  inforaation  is  defined  by  the  row- 
coluan  coordinates  of  the  upper-left  and  lower-right 
corners.  (The  coordinates  are  deteralned  with  an  OCTEK 
function  as  shown  in  Figure  3.2.)  SPIN  transfers  the 
corners  of  the  quartered  teaplate  iaage  Into  the  opposing 
corners  of  a  black  video  file  (see  Figure  3.3).  The  final 
video  file  of  the  teaplate  la  shown  in  Figure  3.4. 
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rifart  3.4  Pinal  Video  nl«  of  TaapUta  •€  P-1*. 


3 . 2  TWO  DIMENSIONAL  FAST  FOURIER  TRANSFORM 

The  correlation  of  the  template  with  the  scene  takes 
place  In  the  frequency  domain  rather  that  in  the  spatial 
domain.  This  Is  accomplished  by  the  multiplication  of 
Fourier  transforms.  A  correlation  In  the  spatial  (or  time) 
domain  is  equivalent  to  performing  a  conjugate  multiplica¬ 
tion  in  the  frequency  domeln  (Ref  6).  This  is  shown  in  the 
following  equetion: 

f  (t)*f  (t)  <  — >  F  (f)*F  <f>  (3-1) 

12  12 

where  f  (t)  and  f  (t)  are  time  domain  functions,  F  (f)  and 
12  1 
F  (f)  are  Fourier  transforms  of  those  functions  and  *  is  the 
2 

complex  conjugate  operator. 

Equation  3-1  may  be  written  in  the  two-dimensional 
spatial  domain,  rather  than  in  the  time  domain.  This  would 
yield  equation  3-2: 

f  <i.J>*f  (  1 ,  J)  <  — >  F  <f  >*F  <f>  (3-2) 

1  2  12 

where  now  f  (i,J)  and  f  (l,j)  are  spatial  functions  with  the 

1  2 

coordinates  (i,j).  To  do  this,  the  video  file  must  first,  be 

in  couples  fora.  The  program  VDTOCF  (Video-To-Coaplex)  does 
this  by  creating  a  new  file  which  contains  the  same  data  of 
the  video  file  but  in  complex  fora  with  the  imaginary  parts 
set  equal  to  sere.  In  this  new  format  the  Two  Dimensional 
Fast  Fourier  Transform  (2DFFT)  may  be  performed. 


■*.  i 


. F* 
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It  is 


Ths  2DFFT  is  psrfersisd  by  the  program  DIRECT,, 
this  procedure  which  requires  the  such  computation  time. 
The  transform  of  a  256x256  complex  file  takes  approximately 
11  minutes  on  the  Eclipse  background  (at  least  22  minutes  on 
the  foreground,  depending  on  the  Job  load).  The  procedure 
of  VDTOCP  and  DIRECT  is  performed  on  both  the  file  con¬ 
taining  the  origin-centered  template  and  the  file  containing 
the  scene  to  be  searched. 


3 • 3  FILE  MULTIPLICATION 

To  correlate  the  two  2DFFTs,  a  point-by-polnt  conjugate 
multiplication  ie  performed.  This  is  done  with  the  program 
VCM  ( Var lable-Con Jugate-Mult ipl lcat ion) .  The  result  is  in 
complex  fora. 

An  alternate  program  was  later  used  to  perform  the  file 
multiplication.  This  was  one  which  also  filtered  the  trans¬ 
form  of  the  template  as  it  multiplied  it  with  the  transform 
of  the  scene.  This  is  discussed  later  in  this  chapter  under 
faction  5.6. 


3.4  INVg 
The 

the  2DPPf 


INVERSE 
about  11 
INVERSE  i 


RSE  TWO  DIMENSIONAL  PAST  FOURIER  TRANSFORM 
output  file  of  VCM  contains  the  multiplication  of 
|T  of  the  template  file  and  of  the  scene  file.  The 
la  to  perform  an  Inverse  2DPFT.  The  program 

Coes  this  computation.  This  function  also  takes 
lntues  on  the  Eclipse  background.  The  output  of 
s  still  in  complex  form. 
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The  resultant  output  file  is  now  converted  to  *  pecked 
vldso  fori  with  tha  progran  CPTOVD  (Cosplex'To-Vldeo).  This 
progres  seals*  tha  coaplex  data  and  convert*  It  Into  integer 
values  fron  0  to  13  so  that  It  way  be  displayed  by  0CT8K. 

3.5  OUTPUT 

The  video  file  containing  the  correlation  Results  nay 
now  be  displayed  on  the  viewing  screen  with  OCTCK  (or  by 
held  copy  with  the  prograa  PIX).  figure  3.5  shows  such  an 
output.  Sharp  white  areas  indicate  possible  correlstione 
with  targets.  The  peaks  nay  also  be  falsh  Correlation  with 
high  energy  noise  in  the  scene.! 

Figure  3.6  is  a  photograph  of  a  display  aonitor  showing 
the  sane  correlation  output  (after  the  grey  levels  were 
reversed).  This  is  Included  ip  order  to  conpensate  for  the 
lost  resolution  of  the  PZX  output  after  reduction. 

Due  to  characteristics  of  this  correlation  in  the 
frequency  donaln  the  peaks  will  occur  at  the  location  of  the 
teaplate  origin  durring  the  line-up  between  the  tenplate  and 
the  target.  In  order  that  the  peaks  correspond  with  the 
actual  location  of  the  targets  the  tenplate  and  the  origin 
oust  exist  at  the  sane  location.  It  is  for  this  reason  that 
the  tenplate  lnforaatlon  was  centered  at  the  origin  of  the 


tenplate  video  file. 
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ri|ur«  3.4  Display  Monitor  Output  of  Figure  3.3. 


3.6  rOUKZKR  riLTERIMO 

As  ssntlonsd  ssrllsr  in  this  chsptsr  t hs  tsnplstss  wars 
originally  titan  froa  tha  aultipla  targat  iaagas  baing 
exaalned.  Hovavtr,  in  an  actual  oparating  systaa  thia  would 
not  ba  practical.  Tha  antlra  procass  was  radona  on  this 


task  scana,  but  a  aara  ganaral  P-16 


An  laaga  of  an  P-16  with  tha  saaa  perspective  was 
digitlsad  using  OCTSK.  This  naw  tsaplata  diffarad  slightly 
in  sisa  and  oriantation  coaparad  t#  tha  targata  in  tha 
scana.  Tha  aost  drastic  diffsranca  was  tha  shading.  Tha 
airplanes  in  tha  original  target  scans  wars  shaded  froa  tha 
side  (Figure  2.2).  This  produced  a  sharp  edge  along  tha 
fuselage  of  tha  aircraft  (Figure  2.7)  after  the  Uac'ga  edge 


detector  was  applied.  Tha  naw  generalised  taaplata  was 
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evenly  shaded  so  that  tha  Wedgs  operator  produced  a  noraal 
outline  of  the  F-16.  (see  Figure  3.7b).  As  a  result,  t'-»e 
new  teaplate  had  a  different  shape  than  the  targets  (shown 
in  Figure  3.7a). 

In  order  to  aatch  a  teaplate  with  non-identical  targets 
a  filtering  technique  ft*  used.  Without  the  filter,  a  non- 
optlaua  output  was  produced.  Figure’  3.6  Is  the  output  of 
the  teaplate  of  Figure  3.7b  correlated  with  the  Wedged  scene 
of  Figure  2.7.  A  photograph  of  the  display  aonltor  output 
is  again  included  ,in  Figure  3.9,  to  show  clarity  of  the 
grey  levels. 

In  theory,  the  filter  is  to  be  applied  to  the  teaplate 
transfora  and  then  the  aultlpllcation  of  teaplate  with  scene 
can  tate  place.  However,  this  would  result  In  an  additional 
one-half  aegabyte  file  being  produced  for  each  filtered 
version  of  the  teaplate.  This  filtered  version  is  refered 
to  as  T '  in  equations  3-3a  and  3-3b.  These  equations  show 
the  theoretical  application  of  the  filter: 


T’»  T*F 


( 3-3a) 


0  •  S*T 1 


( 3-3b) 


where  T  is  the  teaplate  file,  F  is  the  filter,  T*  is  the 
filtered  teap.late  file,  6  is  the  scene  file  and  0  is  the 
output  file.  (The  operator,  *,  represents  coaplex  conjugate 
aultlpllcation.) 


ri|ur«  3.9  Display  Monitor  Output  of  Pi|ur«  3. A. 

Ouo  to  tho  proportion  of  conple*  con|ug*co  ■ultlplict- 
tion,  «i  vo 1 1  oo  oil  Multiplication  oubstltutlon  can  bo 
used.  Tbo  result  of  this  substitution  is  |lvtn  in  squatlon 
3-*: 

0  -  **T»r  (3-4) 

This  oporatlon  is  preforaed  vith  tho  progras  PVCM 
(Pi itered-Varlable-Con juga t e-Mu 1 1 1  pi  teat  ion ) .  Tho  prograa 
PVCM  doos  tho  sano  conjugate  nultiplicatlon  as  does  the 
progran  VCM  (Section  3.3)  but  vith  thn  application  of  a 
variable  slsed  filter. 

Tho  filter  used  was  one  vhich  only  Multiplied  tho  dc 
torn  and  a  specified  nunbsr  of  lover  harnonics.  This  is 
done  bscauso  it  is  the  lover  frequency  harnonics  vhich 
contain  the  aeln  inforaation  of  the  scene,  vhile  the  higher 
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harmonic*  contain  only  minor  details  (Ref  4).  Figure  3.10 
shows  the  output  with  the  application  of  a  filter  of  pixel 
size  7x7  (dc  term  and  tirst  throe  harmonics).  It  should  be 
noted  that  this  was  over  filtered  to  the  point  where  only 
the  high  energy  noise  produced  correlation  peaks.  Figure 
3.11  is  the  output  with  a  more  optimum  filter,  pixel  else 
39x39  (dc  term  and  first  19  harmonics).  The  display  monitor 
output  is  shown  in  Figure  3.12.  This  filter  else  was  judged 
to  be  optimum  for  this  type  of  scene.  FVCM  can  operate  with 
any  user  specified  filter  else  from  7x7  pixels  up  to  233x233 
pixels,  in  increments  of  eight  pixele  (7x7,  13x13,  23x23, 
31x31,  etc.).  The  increment  of  eight  was  needed  due  to  the 
software  techniques  used. 

After  filtering,  the  standard  first  pass  techniques 
were  used.  The  output,  such  as  Figure  3.11',  was  used  to 
determine  the  locations  of  potential  targets  for  the  second 

pass  analysis. 


Figure  3.10  Correlation  Output  of  Target  Scene  (Figure  2.7)  with 
General  F-1S  Template  (Figure  3.7b)  ueing  e  7x7 
Pixel  Filter  (dc  and  firet  3  hereonici) . 

(Over  Filtered) 
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Figure  3.11 


Correlation  Output  of  TarSot  tcano  (Flaura  2  7> 
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•  SECOND  pass  correlation  process 

Th«  basic  procsss  described  In  Chapter  Three  la 
repeated  a  second  tine  to  distinguish  the  correlation  peaks 
as  being  targets  or  noise.  Due  to  the  choracterist lea  of 
the  correlation  process,  high  energy  noise  will  create  a 
peak.  The  second  pass  is  used  to  re-exaalne  the  potential 
target  at  the  location  of  the  correlation  peaks..  The  unique 
feeture  of  the  second  pass  is  that  the  video  file  containing 
this  potential  target  is  energy  noradllsed  to  that  of  the 
taaplate  file. 

A.l  PILE  11CENERATI0M 

Por  reasons  of  efficiency,  portions  of  the  Wedged  scene 
were  extracted  froa  the  256x256  video  file  and  placed  into 
saaller  video  files.  This  was  done  for  two  reasons.  Plret, 
certain  areas  in  the  scene  could  be  close  to  wore  than  one 
peak.  This  would  aean  that  this  ares  could  be  operated  on 
by  the  noraalisation  process  sore  than  ones  and  therefore 
produce  artiflcal  errors. 

The  second  reason  wad  the  tine  factor  involved.  .  Since 
each  file  could  contain  only  one  area  which  was  energy 
noraallsed,  a  file  which  contained  aoltiple  peak  would  have 
to  be  noraallsed  and  racorrslated  for  every  peak.  The 
number  of  peaks  for  potential  targets  could  typically  be  as 
high  as  ten.  Considering  the  fact  that  the  total  process 
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took  about  40  ainutei  for  coaputatlon  of  236x236  files,  thla 
size  vii  prohibitive.  When  saaller  file*  were  uted,  the 
coaputatlon  tine  for  the  prograai  DIRECT  and  INVERSE  waa 
cons iderabl 1 y  shorter. 

The  saall  files  were  created  with  the  progras  SMALL. 
This  prograa  reproduced  the  video  inforaation  froa  one  file 
into  the  saaller  packed  video  file.  The  source  file  was  the 
Wedged  iaage  of  the  scene.  The  location  of  the  potential 
target  was  found  by  using  OCTEK  to  determine  the  row-colunn 
coordinates  of  the  corresponding  correlation  peak.  The  else 
of  the  saaller  video  files  was  liaited  to  128x128,  £4x64, 
and  32x32.  This  was  due  to  the  constraints  on  the  fFT 
progreas  DIRECT  and  INVERSE  for  which  there  was  no  evallible 
source  code.  (Atteapts  were  aade  to  recreate  such  progreas, 
but  the  disk  I/O  tiae  on  the  Eclipse  was  too  lengthly.) 

The  typical  site  of  a  saaller  file  was  64x64.  This 
enabled  the  prograaa  DIRECT  and  INVERSE  to  operate  in  22 
seconds.  With  this  turn  around  tiae,  rcevaluetlom  of 
nuaeroua  potential  target  peaks  could  be  aade.  Saaller 
video  files  of  all  potential  targets  were  aade,  as  well  as  a 
sailltr  video  file  of  the  teaplats.  It  was  the  size  of  the 
teaplate  iaage  which  governed  the  sise  of  ail  these  saaller 
files. 

During  this  second  peas,  the  teaplate  Bust  again  be 
centered  at  the  origin  of  the  file.  This  was  accomplished 
by  the  program  VSPIM  (Variable-Spin),  which,  this  tiae  did 
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4.2  ENERGY  NORMALIZATION 


The  key  to  this  second  pees  ie  the  noreallzatlon  of  the 
energy  In  the  video  files  of  the  potential  targets.  After 
energy  norsaliiation,  recorrelation  should  still  produce  a 
peak  for  a  true  target  but  produce  no  peak  for  high  energy 
noise.  The  energy  in  the  video  iaages  is  defined  ee  being 
the  suaaatlon  of  the  squares  of  the  pixel  values 
(intensity).  This  is  shown  aatheaatically  in  equation  4-1: 

2 

ENERGY  -  2TZ  2  (i.J)  (4-1) 

1  J 

where  1  is  the  intensity  of  the  pixel  at  the  point  (i,j). 

The  noraal ization  routine  used  was  one  which  set  the 
energy  of  the  target  file  equal  to  that  of  the  teaplate 
file.  Since  only  integer  files  were  uied,  the  energies 
could  not  be  set  to  unity.  Therefore,  a  pseudo-unity  was 
used  which  was  equal  to  the  energy  of  the  teaplate  file. 
The  software  to  lapleaent  this  routine  is  contained  in  the 
prograa  NORM. 

This  prograa  creates  a  third  file  which  contains  target 
lnforaatlon  that  aay  be  only  partially  noraallsed.  This 
original  target  file  is  nultlplled  by  a  factor  that  is  the 
square  root  of  the  ratio  of  the  teaplate  energy  to  the 
target  energy.  This  is  given  in  equation  4-2: 
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factor  •  _  / -  (  *  -  2  ) 

V  E« 

where  Ep  la  the  energy  of  the  template  file  and  Eg  la  the 
energy  of  the  target  file. 

Each  new  value  la  then  converted  to  an  integer.  Due  to 
round-off  error  in  the  integer  convertion,  the  proceas  la 
autoaiatlcal ly  repeated  until  no  futher  accuracy  can  be 
achieved. 

NORN  will  not  only  decreaae  high  energy  nolee,  thua 
elininating  a  falae  peak,  but  it  will  elao  enhance  the 
energy  of  a  low  energy,  target  and  increaec  the  correlation 
peak.  Figure  4.3  containa  aone  of  the  video  filea  of  figure 

4.2  after  energy  noraallzation.  Note  that  the  .flret  two 
filea  ahow  lltte  change,  while  the  third  shows  a  significant 
reduction  in  intensity  due  to  the  abundance  of  high  energy 
nolae.  Later,  in  Figure  4.4,  it  will  be  ahown  that  recor¬ 
relation  with  a  true  target  doea  in  fact  aharpen  the  peak 
while  nolae  produce!  a  leaa  defined  peak. 

4.3  KBC0RRE1.ATX0N 

The  next  atepa  of  th<  second  pass  correlation  process 
are  very  aiailar  to  those  of  the  first  pass.  This  tine 
everything  is  redone  but  with  two  differences:  First, 
snaller  files  are  used  (typically  64x64),  ana  second,  a 
conplete  recorrelation  la  done  for  each  potential  target 
area. 
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The  smaller  video  files  ere  declared  cosplex  by  the 
program  VTC  (Vldeo-To-Cosplex)  which  operates  on  variable 
sised  files.  This  is  done  to  the  template  file  as  well  as 
the  target  files.  The  2DFFT  is  again  performed  by  DIRECT. 

Each  of  the  resulting  complex  target  files  are  multi¬ 
plied  with  the  complex  template  file  using  the  program  VCM 
(see  Section  3.3).  The  filtering  multiplier*  FVCM,  has  an 
option  to  filter  64x64  complex  files  with  a  first  seven 
harmonic  filter,  but  it  was  not  used  for  reasons  to  be 
explained  in  the  Results  section  (Section  5.4). 

The  output  of  VCM  under  goes  Inverse  Two  Dimensional 
Fourier  Transformation  with  the  program  INVERSE.  The  result 
is  converted  back  to  packed  video  with  the  program  CTV 
(Complex-To-Video)  which  also  operates  on  variable  sized 
files. 

The  final  output  of  the  recorrelation  can  be  displayed 
on  the  video  screen  with  the  program  VID,  saved  with  OCTEK 
and  printed  with  FIX.  Figure  4.4  shows  such  an  output. 
Due  to  the  nature  of  the  program  CTV,  there  will  always  be 
some  type  of  peak  but  its  location  and  uniqueness,  reveals 
the  level  of  correlation  (see  Chapter  Five).  Visual  inspec¬ 
tion  of  these  video  files  indicate  whether  or  not  a  true 
--target  .was  present  or  if  the  first  pass  correlation  peak  was 


created  by  nigh  energy  noise 
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c)  Correletl 
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Figure  4.4d-e  Second  Pass  Correlation  Output. 

d)  Correlation  Output  iron  Figure  4. 2d, 

e)  Correlation  Output  froa  Figure  4.2e. 
(Mo  reduction  of  PXX  output) 
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v.  evaluation 

Visual  inspection  of  the  second  pass  correlation  was 
confusing  at  Uses.  While  a  sharp  peak  at  Che  center  of  the 
file  wee  clearly  di s t ingul shi ble  as  correlation  with  a  tar¬ 
get,  less  sharp  peaks  becase  uncertain  as  to  their  meaning. 
To  lsprove  the  interpretation  of  the  peaka  a  sore  detailed 
analysis  of  the  nature  of  the  files  was  used  to  produce  the 
guidelines  with  would  be  used  in  the  evaluation  process. 

5 . 1  THREE  DIMENSIONAL  REPRESENTATION 

In  order  to  increase  perception  of  the  information 
contained  in  the  saaller  output  files  of  the  recorrelation, 
a  t hr ee-d leans  tonal  representation  was  sought.  The  Signal 
Processing,  Lab  had  access  to  software  which  produced  a  3-D 
plot  on  either  a  Tektronix  4010  graphics  display  terminal  ->r 
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a  Heathkit  H-19  terminal  with  graphics  capibility.  With  the 
eOtO  an  immediate  printout  was  availible  with  a  Tektronix 
4631  hard  copy  unit. 

The  software  package  was  called  PL0T3D.  It  required 
that  the  Information  source  for  the  plot  be  contained  in  a 
file  of  binary  form.  To  obtain  such  a  form  from  the  stan¬ 
dard  packed  video  format,  a  program  called  BX  was  used  to 
convert  to  binary  fora. 

Figure  3.1  shows  the  3-D  representation  of  an  auto¬ 
correlation  performed  on  figure  4.4a.  This  shows  what  the 
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output  data  should  be  in  tha  bast  peasibla  case  for  this 
type  of  laage.  Figure  S.2  shows  the  plot  of  Figure  4.4b, 
which  is  also  considered  a  good  correlation  peak..  Figure 
5.3  la  the  3-D  plot  of  Figure  4.4e  wbt'.u  is  actually  the 
correlation  of  the  teaplate  file  with  a  file  containing 
noise  and  a  partial  target.  One  should  note  that  this  type 
of  output  contains  a  peak  which  la  off  centered  and  not  as 
proainent  as  the  peaks  in  the  previous  two  diagraae. 


>, 
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3.2  HISTOGRAM  ANALYSIS 

Further  analysis  was  done  on  these  files  by  exaalnlng 
the  histograas  produced  froa  thea.  The  prograa  HISTO 
produced  the  hlstograa  of  a  given  file'.  It  provided  the 
data  in  tabled  fora,  displayed  on  a  terainal,  and  in  graphic 
fora,  displayed  on  a  graphics  tarainal  if  one  were 
availlble. 

Figure  3.4  is  the  hlstograa  produced  froa  the  data  of 
the  auto*correlatl«n  of  Figure  4.4a  (shown  3-D  in  Figure 
3.1).  Since  it  was  an  auto-correlation,  it  was  assuaed  that 
this  1*  the  optiaua  "peaking"  for  this  type  of  laage  data. 
Figure  9.3  is  the  hlstograa  of  Figure  4.4b  (shown  3-D  in 
Figure  3.2).  This  is  the  output  that  can  be  expected  froa  a 
good  correlation  natch  since  the  target  was  auch  like  the 
teaplate  used.  The  hlstograa  of  Figure  9,4  is  froa  Figure 
4 . 4e  (shown  3-D  in  Figure  3.3).  This  is  froa  a  correlation 
which  is  not  judged  to  be  a  good  aatch  of  target  with 
teaplate. 
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mm  fig  4. 4. ft 


Figure  5.4  Hi atograa  of  Pixel  Population  of  Figure  4.4a 


Examination  of  these  three  histograms  indicated  that  a 
correlation  which  was  considered  to  be  good  will  produce  a 
small  number  of  pixels  in  the  high  intensity  regicn; 
whereas,  a  correlation  which  was  judged  to  be  bad  did  not. 
A  bad  match  seemed  to  produce  a  more  even  distribution  of 
pixel  population. 

5.3  EVALUATION  PROCESS 

The  reccrrelated  video  files  were  evaluated  by  con¬ 
sidering  the  two  properties  previously  mentioned.  This  was 
accomplished  with  the  program  EVAL  which  produced  two  scores 
corresponding  to  these  two  properties. 

The  first  score,  Scorel,  rated  distance  of  the  location 

% 

of  the  peak  from  the  center  of  the  file.  A  good  correlation 
produces  a  peak  very  close  to  the  file  center,  due  to  the 
initial  alignment  of  the  first  pass  correlation  with  the 
actual  targets.  Therefore,  the  indication  of  a  good  target- 
template  match  decreases  as  the  distance  to  the  peak  from 
the  center  Increases.  The  peak-to-center  distance  was 
normalized  as  follows: 

D 

N  - - -  (5-1) 

2F 

where  D  is  the  peak-to-center  distance,  F  is  the  file's 
pixel  size  and  N  is  the  normalized  distance.  This  norma¬ 
lized  distance  was  rated  to  produce  Scorel  with  the  function 
in  Figure  3.7. 
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Normalized  Paak-to-C«nt«r  Distance,  N 
Figure  S.7  Function  of  the  Scorel  Evaluation. 

The  aecond  acore,  Score2,  rated  the  eharpneea/unlque- 
neae  of  the  peak.  A  good  correlation  produeea  a  aingle 
aharp  peak.  A  aeaiure  of  the  peak' a  aharpneaa/unlqueneea 
vaa  done  by  a  type  of  integration  of  pixel  population  within 
the  hiatograa.  The  nuaber  of  plxela  in  each  Intenelty  level 
vaa  integrated  froa  level  IS  backvarda  to  zero.  When  5 
percent  of  the  total  nuaber  of  plxela  ^ai  exceeded,  the 
preaent  intenelty  level  waa  recorded.  To  that  level  vac 
added  the  previoua  nuaber  plxela  divided  by  the  percent  of 
the  total  nuaber  which  waa  being  nought. 

Thla  aethod  produced  a  aoaewhat  continuous  output. 
Thla  output  (which  varied  froa  1  for  good,  to  16  for  bad) 
waa  weed  to  produce  Score2.  The  function  which  deteralned 
thla  acore  la  given  in  equation  5-2.  In  thla  equation  the 
output  of  the  integration  deecrlbed  la  refered  to  ae  X.  (X 
la  the  variable  SPOT  in  the  prograa  EVAL.) 
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;  for  X  g  6 

i  for  6  <  XS  16  (5-2) 

The  r'ograe  which  calculated  theae  two  scores,  EVAL, 
aleo  produces  a  third,  total  score  refered  to  as  SCORE. 
Scorel  and  Score2  are  a  rating  froa  0,  for  bad,  to  1 ,  for 
good.  SCORE  is  a  multiplication  of  these  two  factors  (see 
equation  5-3).  Due  to  this  multiplication,  any  rating  <ie- 
ficency  in  either  the  peak  location  or  the  peak  sharpness/ 
uniqueness  will  be  reflected  in  the  total  rating,  SCORE. 

SCORE  -  ( SCORE  1 ) x( SC0RE2 )  (5-3) 

Table  I  is  the  hard  copy  output  of  this  evaluation 
routine.  This  program  also  examines  any  Inconsistency  in 
the  two  scores.  If  the  peak  is  sharp,  but  off  centered,  it 
is  an  indication  of  an  partial  off-centered  target  and  is  so 
stated  in  the  output.  Note  that  EVAL  also  Tanks  the  final 
score  from  one  to  ten. 


TABLE  I 


Evaluation  of  Second  Pat*  Correlation  Output  Piles 


EVALUATION  RESULTS 
SCENE  FILENAME:  FIGURE  4.4 
FILE  NAME  SC0RE1  SC0RE2  TOTAL  SCORE 


FIC404A.  VD 
FI0404B.  VD 
FIG404C.  VD 
FI0404D.  VD 
FI0404E.  VD 


933709 

908891 

889515 

933709 

374610 


1.  000000 
.  971016 
. 961668 
. 956321 
.  666951 


POSSIBLE  OFF-CENTERED  PARTIAL  TARGET 


.933709 
.  882548 
.  855417 
.  892926 
.249847 


RANK 


9 

9 

3 


Note  that  the  first  4  potential  targets  were  rated 
high ,  while  the  fifth  potential  target,  which  was  actually  a 
partial  target  and  high  energy  noise.  Was  rated  such  lower. 
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5.*  RESULTS 


The  process  described  in  this  thesis  itesed  to  work 
quite  well  with  the  lasses  tested.  Appendix  E  contslns 
several  of  the  sdditionsl  iaages  used  for  testing  along  with 
the  results  and  evaluations. 

As  aentioned  earlier  in  Chapter  Pour,  the  aultiplylng 
filter  progrea,  FVCM,  was  able  to  filter  44x64  cited  files. 
However,  it  was  not  used  in  this  capacity*  The  propose  of 
the  filter  was  to  find  specific  targets  with  a  general 
teaplate.  This  was  very  useful  during  the  first  pass  corre¬ 
lation  process  which  purpose  was  to  find  the  location  of 
possible  targets.  It  was  the  function  of  the  second  pass 
correlation  process.  using  energy  noraalisation,  to  re- 
exaalne  the  location  after  the  potential  targets  were  found. 
Trials  showed  that  filtering  during  the  second  pass  seeaed 
to  aerely  blur  the  correlation  output  as  shown  in  Figure 
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Figure  5.S 


Second  Poe*  Correlation  Output  Filee,  With  and 
Without  Filtering. 

a)  Figure  A. 4a,  b)  Figure  4.4a  with  Filtering, 
c)  Figure  4.4e,  d)  Figure  4.4e  with  Filtering. 
(No  reduction  of  FIX  output) 


6.1  CONCLUSIONS 

The  evaluation  process  discussed  in  Chapter  Five  pro¬ 
duced  tangible  results  which  are  Included  in  Appendix  E .  in 
the  cases  where  noise  was  not  present.  or  low  compared  to 
the  energy  of  the  edges  bf  the  target.  the  process  worked 
very  well.  In  trials  where  general  templates  were  used  to 
find  non-identical  specific  targets.  the  filtering  routine 
proved  to  be  effective. 

There  were  exawples  of  the  process  not  being  able  to 
extract  all  of  the  targets  out  of  the  scene.  One  such  case 
was  the  scene  used  throughout  the  text  of  this  thesis. 
Although  the  tewplate  of  the  F-16  extracted  fros  the  scene 
(Figure  3. A)  could  find  all  the  targets  in  the  scene  (Figure 
2.7),  the  general  tewplate  of  Figure  3.6  was  unable  to 
locate  the  two  F-16*s  at  the  outer  edges  of  the  scene. 
(This  case  is  foraelly  excluded  with  Assumption  f 4 ,  Section 
1.4.)  Also,. the  aircraft  just  above  the  srtlfiesl  noise  was 
not  as  clear  as  the  others.  (See  Figure  3*7) 

The  second  pass  correlation  with  energy  norwalistlon 
proved  to  be  extremely  helpful  in  eliminating  the  false 
target  identification  that  occured  due  to  high  energy  noise. 
Sven  in  cases  where  target  correlation  seemed  to  be  somewhat 
difficult,  the  program  SVAL  was  able  to  make  good 
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evaluations  between  the  relative  outputs  given  (see  Figure 
E-l  end  Table  E-II). 

4.2  EECONENDATIONS 

Due  to  tine  constraints,  the  process  developed  in  this 
thesis  has  had  a  finite  asiount  of  testing.  Further  testing 
■ay  be  done  with  higher  background  noise  levels,  less 
■etching  between  target  and  template,  as  well  as  a  combina¬ 
tion  of  these  two  variations. 

Another  possible  advancement  with  this  algorithm  could 
be  its  implementation  with  a  master  program  which  would 
direct  the  functione  of  all  programs  described.  The  visual 
inspection  of  correlation  peaks  after  the  first  pass  could 
be  replaced  with'  software  techniques  similar  to  those  devel¬ 
oped  by  Hamadani  (Ref  2).  Hamadani  developed  a  blob  detec¬ 
tion  technique  which  could  be  modified  to  automat ical ly 
return  the  peak  coordinates  to  a  master  program. 

4. 1  PROCESS  SUMMARY 

The  process  developed  and  described  in  this  thesis  is 
summarised  in  the  flowchsrt  of  Figure  4.1.  This  flowchart 
illustrates  the  eequence  of  every  step  taken  in  this  scene 
anelyels  process. 

The  actual  fils  information  is  described  in  normal 
typo.  The  program  which  produces  the  information  is  in 
capital  letters  and  any  possible  substitutions  are  after 
the  slash.  Tha  aptional  processing  is  in  square  brackets. 
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Figure  4.1a  Flowchart  of  Process  (Pert  I) 
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Figure  4.1b  Flowchart  of  Process  (Part  II) 
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APPENDIX  A 


Material*  and  Eguipaent 

HARDWARE 

Tha  Signal  Processing  Lab  at  the  Air  Force  Institute  of 
Technology  at  Wright-Patterson  AFB,  Ohio  had  all  the  neces¬ 
sary  hardware  and  support  needed  for  this  thesis.  The  I  s.t> 
has  a  Data  General  Nova  2  coaputer  and  a  Data  General 
Eclipse  S/250  coaputer  as  well  as  several  hard  disks  for 
asss  storage.  All  necessary  terainals  and  printers  were 
available. 

The  iaage  Input/output  was  all  done  through  tha  Data 
General  Nova  coaputer.  The  iaage  processing  station  in  the 
lab  contains  an  Octek  2000  Iaage  Analyser,  a  Dage  650  video 
caaera,  an  Electrohoae  video  aor.itor  and  a  Tektronix  A632 
video  hard  copy  unit. 

'The  graphics  station  provided  two  Tektronix  A010 
graphics  terainals,  controlled  by  the  Data  General  Eclipse 
coaputer  and  a  Tektronix  6631  hard  copy  unit  which  was 
connected  to  both  graphics  terainals. 

SOFTWARE 

The  Signal  Procasslng  Lab  had  access  to  several  soft¬ 
ware  packages  which  were  extensivly  used  in  this  thesis. 
The  following  pages  are  explanations  on  the  use  of  tha 
software  for  which  the  source  code  was  unavalble  (Direct, 
Inverse,  and  Plot3D)  and  for  OCTEK  which  was  drove  the  Octek 


Notes  on  the  program  DIRECT 

Use:  Two  Diaemisnal  Discrete  Fourier  Transform  (2DFFT) 

Execution  format: 

DIRECT  inputfile/1  #/N  { ou t put f i 1 e /O J 

Parameter  explanation: 

I ngu t f 1 1 The  input  file  must  be  in  binary  format  since  a 
RDBLK  is  used  to  read  data  from  the  input  file.  The  input 
file  is  assumed  to  contain  complex  data,  with  the  imaginary 
part  sero.  The  input  array  is  assumed  to  be  square,  with  a 
minimum  row  or  column  length  of  64.  If  an  array  of  sisa 
<  64  is  to  be  used  as  input,  each  column  should  be  stored  as 
one  block  of  data  even  though  only  the  firft  N  complex 
numbers  in  the  block  correspond  to  actual  data  points.  Only 
N  rows  need  be  stored.  Thus,  for  N  <  64,  ths  size  of  the 
input  array  must  be  64  by  N  columns. 

#2  An  Integer  must  be  specified  to  indicate  the  size  of 

the  input  ,  array  -(NxN),  which  is  ths  ilia  of  the  DPT 
computed.  N  must  be  a  power  of  two,  N  <  512. 

Outgutfile^  The  output  file  specification  is  optional.  If 
no  output  file  is  specified,  the  transformed  data  will  be 
written  back  into  the  input  file,  destroying  the  input  data. 
The  output  file  is  complex  and  in  binary  format.  The  origin 
of  the  DPT  samples  will  be  at  the  location  (N/2,  N/2),  and 

the  rows  and  columns  will  be  transposed.  Note:  This  pro¬ 

gram  runs  much  faster  if  the  outputfile  is  not  specified. 
(Ref  6) 
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Note*  on  the  program  INVERSE 

Use:  Two  Dimensional  Inverse  Discrete  Fourier  Transform 
Execution  format: 

INVERSE  input file/I  ( out put f i 1 e /0 )  «/N 
Parameter  explanation: 

Ingutfile^  The  input  file  must  be  in  binary  format  aincc  a 
RDBLK  is  used  to  read  data  from  the  input  file.  The  input 
file  is  assumed  to  contain  complex  data,  and  must  have 
minimum  length  of  64.  If  the  NxN  array  has  N  <  64,  the 
columns  should  be  augmented  with  zeros  so  that  each  column 
can  be  stored  as  one  block  of  data.  Thus,  for  N  <  64,  the 
size  of  the  input  array  must  be  64  rows  by  N  columns  of 
complex  DFT  samples.  The  origin  of  the  frequency  axes  must 
be  at  (N/2  ♦  1,  N/2  +  1). 

Outgutfile^  The  output  file  specification  is  optional.  If 
no  -output  file  is  specified,  the  inverse-transformed  data 
will  be  -written  back  into  the  input  file,  destroying  the 
input  data.  The  output  file  is  complex  and  in  binary  for¬ 
mat.  Note:  This  program  runs  much  faster  if  the  outputflle 
is  not  specified. 

#2  An  Integer  must  be  specified  to  indicate  the  size  of 
the  input  array  (NxN),  which  is  the  size  of  the  DFT  com¬ 
puted.  N  must  be  a  power  of  two,  N  <  312.  (Ref  6) 
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Notes  on  pregrts  PL0T3I 

Use:  Three-Dimensional  Representation  of  a  Binary  Pile 

Execution  format: 

/ 

PL0T3D  input  f 1 le/ I 

Explalnation: 

Ingu t f 1 1 e^  The  input  file  to  this  program  must  be  written 
in  binary  format.  This  is  accomplished  on  video  files  with 
the  program  BI. 

This  program  must  be  run  on  a  terminal  with  graphics  capa¬ 
bility,  such  as  a  Tektronix  4010  or  an  H-19  with  a  graphics 
package .  A  hard  copy  ia  avalllble  via  a  Tektronix  4631 
tardcopler.  The  inputs  prompts  are  self-explanitory  with 
:he  exception  of  the  following: 

XSIZE:  Relative  base  size  in  the  X  direction  on  the 

screen  (typically  8). 

YSIZE:  Relative  base  size  in  the  Y  direction  on  the 

screen  (typically  3). 

HEIGHT:  Relative  size  in  the  Z  direction  (typically  3). 

NORMALIZATION:  Yes,  the  output  of  BI  needs  to  be  norma- 


iota*  on  th«  pro|ria  OCTEK 

Ise:  Operation  of  the  Octek  2000  Iaage  Analyzer 

Execution  foraat : 

OCTEK 

xplanation: 

he  following  Menu  le  the  available  options  in  the  prograa 
CTEKs 


tSV:! 


1 . 

ACQUIRE 

'  - 

2. 

L0A0 

- 

display  video  (.VD)  file 

3. 

SAVE 

- 

aave  diaplay  in  viaeo  (.VD)  file 

A. 

CURSOR 

- 

activate 

cureor  ot  display  pixel 

values 

5. 

CLEAR 

- 

clears  display 

6. 

BARS 

- 

display  greyscale  bar  test  pattern 

1 1  . 

ACQUIRE2 

- 

take  picture  (with  TV  caaera  #2) 

1  2. 

BINARIZE 

- 

threshold 

display 

13. 

LOOKUP 

- 

liter  out 

put  translation  table 

1A. 

NECATIVE 

- 

display  n 

egatlve 

ii. 

GREYSCALE 

- 

restore  t 

rue  greyscale 

r 

16. 

FILL 

- 

fill  bloc 

k 

-1  . 

QUIT 

- 

exit  prog 

TEK  optione  1,  2  and  3  wer 

1  video  files.  Option  A, 


e  used  for  input  and  output  of 
CURSOR,  was  used  to  detersine 


e  location  of  -correlation  peaks.  Option  12,  BINARIZE,  was 
ed  to  anaylise  correlation  output  files.  Option  16,  FILL, 
e  used  to  eliainate  unwanted  noise  in  a  teaplnte  video 
le. 
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Note*  on  libary  PICBUF 

Use:  Virtual  Neaory  Syitea  for  Digital  Picture  Processing 

Reloader  reference: 

FSPICBUF.LB  (for  Fortran  V  on  the  Eclipse) 
KF4PICBUF.LB  (for  Fortran  IV  on  the  Nova) 

Explanation: 

PICBUF  is  a  libary  of  Fortran  subroutines  which  implements  a 
virtual  memory  system  capable  of  storing  and  accessing  two- 
dimentional  digits!  data.  The  PICBUF  routines  crests 
virtual  memory  buffers,  acceeee  image  date  rows  end  pixel 
points,  and  store  the  processed  buffers  into  disk  files.  An 
extended  explanation  of  these  routines  is  Included  in 
Beference  5. 

The  PICBUF  libary  was  used  in  the  following  programs  and 
subroutines: 

SPIN.PE 
MOVE. Ft 
BMALL.FR 
VBPIN.FR 

vid.fr 

windov.fr 

■val.fr 
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Notes  on  prograa  PIX 


Use:  Hardcopy  output  of  video  file. 

Execution  foract: 

PIX 


Explanation: 

This  pro|r«e  will  convert  video  pixels  to  lineprlnter 
pixels,  and  will  put  the  picture  in  a  file  or  to  the 
Prlntronlx  300  lineprlnter.  This  pro|raa  prints  either  a 
coaplete  236x236  pixel  picture,  or  a  saeller  picture  that  is 
of  a  desired  length  specified  by  the  user. 


This  prograa  was  originally  written  by  Lt.  Sisaons,  revised 
by  Lt.  Croaer  and  later  by  Jla  Vets.  (Ref  1,9) 
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APPENDIX  B 


Preprocessing  Software 

This  Appendix  contains  the  prepocesslng  software  used 
to  perform  edge  extraction  on  the  video  files.  Included  are 
the  following  programs  and  subroutines: 

WEDCE.FR 

WR0W1.FR 


COMB.FR 


o  o  u 


£**•**•**•**•**•*«*•*****•***«*»•* ft***********.**************** 

c 

C  WEDCE.FR  -  DC  FORTRAN  5  -  CAPT  ROBERT  WELLS,  DEC  1983 
C  REDONE  BY  LT  RICHARD  MILLS  JULY  1984 

C 

C  This  program  edges  the  input  file  (packed  video  format) 

C  with  a  mask  operator  called  the  Wedge  Operator.  The 

C  output  is  separated  into  four  256  X  256  packed  video 

C  files  according  to  which  mask  produced  the  maximum  value. 

C  If  two  masks  of  different  orientations  produce  the 

C  maximum,  that  value  goes  to  both  output  files. 

C 

C  RELOAD  COMMAND: 

C  RLDR  WEDCE  UNPACK  REPACK  WR0W1  ?FLIB@ 

C 

C  EXECUTION  COMMAND: 

C  WEDCE 

C 

c*«s* seseeseeeseeeeseesaes* **********  ******* *********  ********** 


INTECER  I NMAT (256, 1 6 > . TEMP ( 5 1 2 )  ,  TOCCLE 

I NTECER  OMAT 1(256.8), 0MAT2 (256,8), 0MAT3 ( 256,8). OMAT4 (256.8) 
INTECER  INFLNM( 7 ) , OFLNM 1 (7) ,OFLNM2(7) , 0FLNM3 ( 7 ) »  0FLNM4 ( 7 ) 

C 

C  ACCEPT  INPUT 
C 

ACCEPT"ENTER  INPUT  FILENAME  ->"  4? 

READ ( 1 1  , 1)1 NFLNM( 1  ) 

1  FORMAT ( S 1 3 ) 

ACCEPT "ENTER  OUTPUT  FILENAME  #1  ->" 

READ ( 11,1 )OFLNM 1 ( 1  ) 


ACCEPT" ENTER  OUTPUT  FILENAME  #2  ->" 
READ( 11,1 ) 0FLNM2 ( 1 ) 

ACCEPT" ENTER  OUTPUT  FILENAME  #3  ->" 
READ( 1 1 , 1 )OFLNM3( 1 ) 

ACCEPT "ENTER  OUTPUT  FILENAME  #4  ->" 
READ ( 11,1 )OFLNM4( 1 ) 


CALL  OPEN( 1 , INFLNM, 1 , IER) 

CALL  CHECK(IER) 

DELETE  OUTPUT  FILES  IF  THEY  EXIST  AND  (RE)CREATE  THEM 

CALL  DFILW  (OFLNM 1,1 ER ) 

CALL  CFILW  ( OFLNM 1 ,2 , KER ) 

CALL  CHECK ( KER  ) 

OPEN  2,  OFLNM 1 , ATT • "OR" , ERR* 100 
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CALL  DFILV  ( 0FLMM2  ,  I  ER  ) 

CALL  CFZLW  (0FLMM2 , 2 , KIR) 

CALL  CHECK(KER) 

OPEN  3.  OFLNN2  ,  ATT* "OR" , ERR* 100 

CALL  OFILW  (0FLNM3 , IER) 

CALL  CFZLW  (0FLNM3 , 2 , KER ) 

CALL  CHECK(KER) 

OPEN  A,  0FLNM3 , ATT* "OR" , ERR* 100 

CALL  OFZLV  (0FLNM4.IER) 

CALL  CFZLW  ( 0FLNM4 , 2 , KER ) 

CALL  CHECK(KER) 

OPEN  3,  OFLMM4,ATT-"OR",ERR-100 

DO  2  1*1,25*  i  FZLL  FZRST  AMD  LAST  4  ROWS  OF 
TEMP(I)-0  ;  OF  OUTPUT  PZLES  WITH  ZERO 
2  COMTZMUR 

CALL  WRRLK( 2.0, TEMP, 1 , ZIR) 

CALL  CHRCK(ZER) 

CALL  VRBLK( 2 , 43 , TRMP , I , ZER ) 

CALL  CHECK(IER) 

CALL  WRILK(3. 0, TEMP, 1 , ZER) 

CALL  CHECK(IEK) 

CALL  VRILK( 3,63 .TEMP , 1 , Z ER) 

CALL  CHECK(ZER) 

CALL  WRRLK(4, 0, TEMP, 1 , IER) 

CALL  CHECK(IER) 

CALL  VRBLKC  4,63, TEMP , 1 , I ER ) 

CALL  CHECK(IER) 

CALL  WRRLK( 3 .O.TEMP, 1 , ZER) 

CALL  CHECK(IER) 

CALL  WRBLK( 3,43 .TEMP , 1 , ZER) 

CALL  CHECK(IER) 

CALL  RDBLK<I,0,TEMP,2,IER)  {READ  ZH  FIRST  •  ROWS 
CALL  CHECK(IER) 

CALL  UNPACK(312 .TEMP, ZNMAT) 

TOGGLE-O  ;  FLAC  TO  SHOW  BUFFER  WRAP-AROUND 

DO  3  1*1,31  ;  COUNTER  FOR  2K  BUFFER  LOADS 

WRITE  FREE ( 10)1  ;Kaap  track  of  running  program 

CALL  RDBLK( 1 , 2*1 .TEMP, 2, IER)  ; READ  NEXT  S  ROWS 
CALL  CHECK(IER) 

CALL  UNPACK (312, TEMP , I NMAT ( 1 , 9-B*T0CGLB ) ) 

DO  3  K* 1 , S  ;  COUNTER  FOR  OUTPUT  BUFFER  ROWS 

CALL  WR0W1 (K, TOGGLE, INMAT.0MAT1 , 0MAT2 , 0MAT3 , 0MAT4 ) 
CONTINUE 
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CALL  REPACK( 5 1 2 .OMATl .TEMP) 
CALL  WRBLK (2, 2*1-1 ,TEMP,2,IER) 
CALL  CHECK ( I ER ) 

CALL  REPACK ( S 1 « , OMAT2 , TEMP ) 
CALL  WRBLK(3, 2*1-1 , TEMP, 2.IER) 
CALL  CHECK ( I ER  ) 

CALL  REPACK ( 3 1 2 , OMAT3 , TEMP ) 
CALL  WRBLK,(  4, 2*1-1  ,  TEMP  ,2,1  ER  ) 
CALL  CHECK ( I ER ) 

CALL  REPACK ( 512, OMAT4 , TEMP ) 
CALL  WRBLK<3, 2*1-1 .TEMP, 2, IER) 
CALL  CHECK (IER) 


WRITE  8  ROWS  OF  RESULTS 


IF(TOCGLE.EQ. 1 ) GO  TO  4 

TOGGLE- 1  ;  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAG 

CO  TO  V 

A  TOGGLE-O 

5  CONTINUE 


CALL  RESET 

STOPh<7><7><7X7>WEDCEm 

END 
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SUBROUTINE  VROW1 ( K , TOGGLE , I NMAT, OMAT 1 , OMAT2 , OMAT3 . OMAT4 ) 

BY  WELLS  1983 

INTECER  K, TOGGLE, INMAT(256, 16) .HEDGE 

I NTECEK  OMAT 1(236,8). OMAT2 ( 236,8), OMAT3 (236,8), OMAT A ( 236,8) 
REAL  A1 ,A2,A3,A4,B1 ,B2,B3,B4,C1 .C2.C3.C4.D1 ,D2 , D3 , D4 , MAXIM 

I F ( TOGGLE .EQ. 1 ) GO  TO  1  ;  FIND  INPUT  BUFFER  INDEXES 

Kl«K+3 
K2-K+4 
K3-K+3 
GO  TO  2 

1  K1-MOD(K+10, 16)4-1 
K2»MOD( K+l  1,16)4-1 
K3-NOD(K+12, 16)  +  1 

2  CONTINUE 

DO  3  J-1,4  ;  ZERO  FILL  FIRST/LAST  4  COLUMNS  OF  RESULTS 

OMAT1(J,K)-0 
0MATl(J+232,K)-0 
0MAT2(J,K)-0 
OMAT2( J+252 ,K)-0 
OMAT3( J,K)-0 
OMAT3( J+252 , K)“0 
OMAT4(J,K)-0 
OMAT4( J+252 ,K)*0 

3  CONTINUE 

DO  4  J-5,232  ;  EDGE  THESE  COLUMNS 

Jl-J-1 
J2*  J 
J3-J+1 

Al»0.600*( INMAT ( J 1 , K 1)+INMAT( J2 , K 1 )  +  INMAT( J3 ,  K 1  )  )  - 
+  0.3OO*(INMAf( J! ,K2)+INMAT(J2,K2>+XNMAT(J3,K2>+ 

+  INMAT (  J 1  ,  K3  )  +  INMAT(  J  2  ,  K3  )  + INMAT (  J3  ,  K3  )') 

A 2  *0 . 600* ( INMAT (J1 ,K3)  +  INMAT(J2,K3 )  +  INMAT( J3 ,  K3 ) ) - 
+  0.300*(INMAT(J1 ,K2 )+I NMAT (J2,K2)+INMAT(J3,K2)+ 

♦  INMAT( J 1 , K 1 )+INNAT( J2 , K1 )+XNMAT( J3 , K1 ) ) 

A3»-A  1 
A4--A2 

Bl»0.600*( INMAT ( J2 , K 1 )+INMAT( J3 , K 1 )+INMAT( J3 ,K2 ) )- 

♦  0.300*(INMAT(J1,K1 )+INMAT(J2,K2)+INMAT(J3,K3)+ 

+  INMAT ( J 1 , K2 )+I NMAT ( J 1 , K3 )+INMAT (J2,K3) ) 
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B2-0. 600*(INMAT< J1 , K2 )+INMAT(Jl . K3 )+INMAT(J2 ,K3> >- 
0. 300* ( XNMAT( J1 , K1 )+INMAT( J2 , K2 ) +INMAT( J3 , K3 )+ 

I NMAT ( J  2 ,KI  )  +  INMAT(J3 , K1  )  +  I NMAT( J3 , K2 ) ) 

B  3  *  -  B  1 
B4--R2 

Cl«0.600*< I KM AT ( J3.K1 )+I NMAT ( J 3 , K2 )+I NMAT< J3 , K3 ) ) - 
0 . 300*  (  I  NMAT  (  J2,K1)+INMAT(J2  ,  K2  >+INMAT( J2 , K3 )♦ 

I NMAT ( J 1 ,K1 ) + INMAT ( J I , K2 )+I NMAT< J1 , K3 ) > 

C2»0.600*( I NMAT ( J 1 , K 1 )+I NMAT ( J 1 , K2 )+INMAT< J1 ,K3))- 
0 . 300*  (  INMAT  (J2,K1)  +  X  NMAT  (  J2.,  K2  )+INMAT( J2  ,  K3  )+ 
INMAT( J3,Kl )♦ INMAT ( J3 , K2 )+I NMAT( J3 , K3 ) ) 


C3--C1 

CA»-C2 

Dl-0. 600*  <  INMATCJl  ,K2  )+lNMAT( J1  ,  K  1  )♦  I  NMAT<  J2  ,  K 1  ))- 
0 . 300* ( I NMAT ( J 1 ,K3)  + I NMAT ( J2 , K2 )+INMAT(J3 , K1 )+ 
XNHAT( J2 ,K3)+INMAT< J3 ,  K3)  +  INMAT( J3,K2)> 

D2-0.600*(INMAT(J2,K3)+INMAT(J3,K3)+INMAT(J3,K2))- 
0. 300*<INMAT(J1  ,K3)+INMAT( J2 ,K2)+INMAT(J3,K1  )■*■ 

I  NMAT  (  J  1  ,  K2  )+ INMAT (J  1  ,K1  )<f  XNMAT(  J2  , kl  )  ) 

03  — D1 


DA  — D2 

MAXIM-AMAX1 ( A  1 , A 2 , A 3 , A A , B 1  , B2 , B3 , B4 , C 1 , C2 , C3 , C4 , D 1 , D2 , D3 
,04) 

OMAT1(J,K)»0 
OMAT2(J,K)-0 
OMAT3 ( J , K ) “0 
0MAT4(J,K)-0 

MfcDGE«»ANINT(  MAXIM)  ;  INTEGRAL  EDGE  VALUE 

IF(MEDCE.CT. 1 5)MEDCE-l 5  ;  TEST  FOR  OVERFLOW 
IF(MEDG£.LT.0)MEDCE-0  ;  AND  UNDERFLOW 

IF(Al . EQ. MAXIM. OR . A2 . EQ . MAX XM . OR . A3 . EQ . MAXIM .OR . A4 . EQ . 
MAXIM)0MAT1 ( J,K)-MEDGE  ;  ORIENTATION  fl 

I F( B 1 . EQ. MAXIM. OR. B2-E0. MAXIM. OR. B3 . EQ . MAXIM . OR . B4 . EQ . 
MAXIM ) 0MAT2 (  J  ,  K  ) ■HEDGE  ;  ORIENTATION  #2 

I F ( C 1 . EQ. MAXIM. OR. C2 . EQ . MAX IM . OR . C3 . EQ . MAXIM. OR . C4 . EQ. 
MAX IM)0MAT3 ( J , X ) “MEDCE  ;  ORIENTATION  f3 
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I F  <  0 1 . EQ . KAX IM.0R.D2.EQ. MAX I M . OR . D  3 . EQ . MAX I M . OR . D4 . EQ . 
♦  MAXIM )OMAT  4  (  J ,K)"MEDGE  ;  ORIENTATION  #4 

4  CONTINUE 

RETURN 

END 


C  PROGRAM  COMB  - 

DC  FORTRAN  5  -  BY  CAPT  ROBERT  WELLS,  DEC  19B3 

C  THIS  PROGRAM  COMBINES  UP  TO  FOUR  VIDEO  IMAGES  (PACKED  VIDEO 
C  FORMAT).  EACH  PIXEL  OF  THE  INPUT  IMAGES  IS  TESTED  TO  SEE  IF 

C  IT  IS  NON-ZERO .  THE  VALUE  OF  THE  FIRST  NON-ZERO  PIXEL  OF  THE 

C  INPUT  IMAGES  IS  ASSIGNED  TO  THE  OUTPUT  IMAGE  AT  THAT  PIXEL 

C  LOCATION.  IF  ALL  INPUT  IMAGES  ARE  ZERO,  THE  OUTPUT  IS  ZERO. 

C  THE  OUTPUT  IS  STORED  IN  PACKED  VIDEO  FORMAT. 

INTECER  INFNM1 (7) , INFNM2 ( 7 ) , I NFNM3 ( 7 ) , INFNM4(7) ,OUTFNM(7) 
INTEGER  INI  (2048) ,SN2< 2048) ,IN3(204B) , IN4( 2048) .OUT ( 2048) 
INTEGER  TEMP( 512), FLCNT 

1  FORMAT ( S 1 3 ) 

2  ACCEPT "ENTER  NUMBER  OF  FILES  TO  BE  COMBINED  ->", FLCNT 
IF< (FLCNT.LT. 2) .OR. (FLCNT. GT. 4) )GO  TO  2 

CO  TO  (2,4,3)  FLCNT 

ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ ( 1 1 , 1 )INFNM4( 1 ) 

CALL  OPEN( 5 . INFNM4 , I , IER) 

CALL  CHECK(IER) 

3  ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ( 11,1) INFNM3 ( 1 ) 

CALL  OPEN(4,INFNM3,i,IER) 

CALL  CHECK(IER) 

4  ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ(ll.l) INFNM2 ( I ) 

CALL  OPEN( 3 , INFNM2 , i , IER) 

CALL  CHECK(IER) 

ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ ( I  1 , 1 ) INFNM1 ( 1 ) 

CALL  OPEN( 2 , INFNMI , I , IER ) 

CALL  CHECK(IER) 

.  «««»«  UUiPUT.  FILENAME  ->" 

KfcAU( 11,1 )OUTFNM( 1 ) 

CALL  OPEN( 1 .OUTFNM, 3, IER) 

CALL  CHECK ( I ER ) 

DO  14  1-0,31  ;  COUNT  FOR  2K  BUFFER  LOADS 

GO  TO  (2,6,5)  FLCNT 

CALL  RDBLK( 5,2*1 , TEMP , 2 , IER ) 

CALL  CHECK(IER) 

CALL  UNPACK(S12,TEMP, IN4) 


SI 


5  CALL  RDSLK(4,2*I,TEMP,2,XER) 

CALL  CHECK ( I ER ) 

CALL  UNPACK( 512 , TEMP , IN3  ) 

6  CALL  RDBLK(3,2*I,T£MP,2,IER) 

CALL  CHECK ( I ER ) 

CALL  UNPACK(512,TEMP,IN2) 

CALL  RDBLK(2,2*I,T£MP,2,XER) 

CALL  CHECK(IER) 

CALL  UNPACK ( 512, TEMP, INI) 

DO  13  J-1,2048  ;  COMBINING  LOOP 

CO  TO  (2,6,7)  FLCNT 

IF(IN4(J).NE.O ) GO  TO  9 

IF(IN3(J).NE.O)GO  TO  10 

IF(IN2( J) .NE.O)CO  TO  11 

IF(IN1 ( J) .NE.O)GO  TO  12 

OUT (  J  )  »0  {ZERO  IF  ALL  IN-FILES  ARE  ZERO 
GO  TO  13 

0UT(J)-IN4< J) 

GO  TO  13 

OUT (J)*IN3(J) 

GO  TO  13 

OUT< J)«IN2< J) 

GO  TO  13 

12  0UT(J)-IN1(J) 

13  CONTINUE 

CALL  REPACK(312 , OUT, TEMP) 

CALL  WRBLK ( 1 ,2*1 .TEMP, 2 , IER) 

CALL  CHECK ( IER ) 

14  CONTINUE 
CALL  RESET 

STOP"<7><7><7><7>COMBH 

END 


7 
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APPENDIX  C 

Correlation^  Normal Izat lon1  and  Evalutlon  Software 

This  Appendix  contains  the  main  software  developed  and 
used  in  this  thesis.  Included  are  the  following  prograsis 
and  subroutines: 

SPIN . FR 

MOVE . FR 

VDTOCP.FR 

VCM.FR 

FVCM.FR 

CPTOVD.FR 

SMALL.FR 

VSPIN.FR 

NORM . FR 

VTC.FR 

CTV.FR 

EVAL.FR 
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c*«** *********** ****** *******  ******** A  ************************* 

c 

C  SPIN.FR  -DC  FORTRAN  5-  LT  RICHARD  L.  MIL  '  27  JULY  1984 

C 

C  This  program  will  accept  the  name  of  a  file  which  contains 

C  a  template  and  create  a  new  file  which  place*  the  template 

C  center  at  the  orgln  of  the  video  file.  The  only  other  info 
C  needed  li  the  X-Y  coordinates  (Row-Column)  of  the  original 
C  template's  position.  The  background  of  the  new  file  will 
C  be  set  to  black  (pixel  value  ■  0). 

C 

C  RELOAD  LINE: 

C  RLDR  SPIN  MOVE  F5PICBUF.LB  @FLIB<* 

C 

C  COMMAND  LINE: 

C  SPIN 

C 

C*********** **********************************  *****************! 

C*  PICBUF  STUFF 

PARAMETER  NBUFSZ-300 

PARAMETER  NCOLSMAX-236 

INTEGER  IBUF(NBUFSZ) , I BUF2 ( NBUFS2 ) 

INTEGER  I ARRAY ( NCOLSMAX ) t I NAME ( 20 ) , XOTNM( 20 ) 

IMPLICIT  INTEGER(A,B,C,R,X,Y) 

C*  I/O  CONSTANTS 

PARAMETER  IN- 11 
PARAMETER  OUT- 10 
C 

C  ACCEPT  INPUT 
C 

ACCEPT "NAME  OF  THE  INPUT  FILE:-->" 

READ ( IN , 40 ) INANE ( 1 ) 

40  FORMAT ( S40) 

ACCEPT"NAME  OF  THE  OUTPUT  FILE: — >" 

READ( IN, 41 )I0TNM( 1 ) 

41  FORMAT(S40) 

C 

C  PICBUFF  FOR  INPUT  FILE 
C 

CALL  PICFMTOBUF,  INANE,  IHDR) 

CALL  GFMT( I BUF  , NR0WS , NC0LS , NB ITS , I MODE ) 

CALL  MAKB(IBUF) 

CALL  PICIN(IBUF, I NAME , IHDR) 

CALL  DFILW  ( I0TNM , I ER )  (delete  output  file  if  it  exists 
1F(IER.NE. 1 )G0T0  90 
90  CONTINUE 

C 

C  PICBUFF  FOR  OUTPUT  FILE 
C 


.V.1 


$3 


i-.v; 


o  o  n  n  <■>  o 


CALL  CFMT< IBUF.NROWS , NCOLS , NBITS , IKODE' 
CALL  PFMT ( I BUF 2 , NROWS .NCOLS, NBITS, I MODE ) 
CALL  MAKB( IBUF2 > 

C 

C  GET  POSITION 
C 

A  TYPE"  " 

TYPE "ROW, COLUMN  OF  THE  UPPER  LEFT  CORNER" 
ACCEPT"Ot  THE  ORIGINAL  TEMPLATE-- >", R , C 
XL-C 
YL«R 


ASSUME  X  INCREASES  TO  THE  RIGHT  AND  Y  INCREASES  DOWNWARD 
TYPE"  " 

TYPE"ROW, COLUMN  OF  THE  LOWER  RIGHT  CORNER" 

ACCEPT"OF  THE  ORIGINAL  TEMPLATE — >",R,C 

XH-C 

YH“R 


VERIFY  PROPER  ORIENTATION 

IF( (XL. LT.XH) .AND. (YL. LT. YH) )  GO  TO  5  ;  OKAY 

TYPE"  " 

TYPE "ORIENTATION  IS  INCORRECT  ...  RE-ENTER - >"  ;  MOT  OKAY 

GO  TO  4 


:  CALCULATE  MID  PT. 

* 

5  XM*»XL+(XH-XL  )  /  2 

YM«YL+(YH-YL>/2 

XM-INT(XM)  ; INSURE  INTEGERS 

YM-INT(YM)  ; INSURE  INTEGERS 


MOVE  UPPER  LEFT  CORNER 

AXL»XL 

AYL-YL 

AXH-XM-1 

AYH-YM-l 

BXH-2S6 

BXL-256-( AXH-AXL) 

BYL»2  56- ( AYH-AYL) 

CALL  MOVE( AXL» AXH , AYL  »AYH ,  BXL , BXH , BYL , I BUF , IBUF2 , 1 ARRAY ) 


MOVE  LOVER  LEFT  CORNER 


AXL-XL 
AXH-XM- I 
AYL-YM 
AYH-YH 

BYL-  * 

BXH-2S6 

BXL«256-(BXH-BXL) 

CALL  MOVE (AXL , AXH  ,  AYL  v  \YH , BXL , BXH , BYL , I BUF , I BUF2, 1 ARRAY ) 


MOVE  UPPER  RIGHT  CORNER 

AXL-XM 

AYL-YL 

AXH*XH 

AYH-YM-1 

BXL-1 

BXH- l+( AXH-AXL) 

BYL-256-(AYH-AYL) 

CALL  MOVE( AXL , AXH  ,aYL(AYH,  BXL , BXH  ,  BYL, I BUF , I  BUF 2 , I ARRAY ) 


MOVE  LOVER  RIGHT  CORNER 

AXL-XM 

AYL-YM 

AXH-XH 

AYH-YH 

BXL-  1 

■  BYL-I  ' _ 

BXH- !♦< AXH- AXL) 

CALL  MOVE(AXL(AXH,AYL.AYH,BXL,Br.H,BYL,IBUF,IBL,F2 , I ARRAY) 

CALL  PICOUT( I BUF 2 , IOTNM, IHDR)  to  flic 

CALL  R2LB(IB1’F)  jclose  buffers 

CALL  2ELB( i BUF 2 ) 

CALL  RESET 

STOP"<7><7X?><7>*PIN" 


c**« ***** ************************ ***************************** 

;  >;  c 

i  '  C  NOVE.FR  -DC  FORTRAN  5-  LT  RICHARD  L.  MILLS  27  JULY  1984 

I  C 

■  C  CALLED  BY: 

;  C  SPIN. FR 

C 

.  c****************** ft****************************************** 

1C  Subroutine  MOVE  will  do  the  actual  pixel  duplication  In 

C  the  output  file  (region  B)  from  the  input  file  (region  A). 
c***a **************************** ***************************** 


I 


I  @ 


SUBROUTINE  MOVE ( AXL , AXH , AYL , AYH , BXL, BXH , »YL , 11UF , IBUF2 , 
+  I ARRAY ) 


IMPLICIT  INTEGER  (A,B,X,Y) 

TYPE"MOVINC  PIXELS" 

ILEN-AXH-AXL+1 

IR0W2-BYL 

ICOL-AXL 

IC0L2-BXL 

DO  200  I ROW* AYL , AYH 

CALL  GROVdBUF,  IROW,  I  COL  ,  I LEN  ,  I  A, '.RAY  > 
CALL  PR0V(IBUF2 . IROW2 , IC0L2 , ILEN . I ARRAY ) 
IROV2-IROW2+1 
200  CONTINUE 
RETURN 
END 


t 
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C  PROGRAM  VDTOCP  -  DG  FORTRAN  3  -  CAPT  ROBERT  WELLS,  DEC  19S3 

C  THIS  PROGRAM  CONVERTS  A  256  X  256  PACKED  VIDEO  FILE  TO  A 

C  256  X  256  COMPLEX  FILE 

INTEGER  IFLNM(7) , OFLNM( 7  ) 

INTEGER  MAIN (7),F3(7),MS(2),S1(2),S2(2),S3(2) 

INTEGER  TEMP (512) ,1N(2046) 

COMPLEX  OUT ( 2048 ) 

CALL  IOF( 2 , MAIN , I FLNM , OFLNM , F3 , MS ,S1,S2,S3) 

CALL  OPEN( 1 , I FLNM, 1 , IER) 

CALL  CHECK ( IER ) 

CALL  OPEN (2, OFLNM, 3, IER) 

CALL  CHECK(XER) 

DO  3  1-0,31 

CALL  RDBLK( 1 ,2*1, TEMP , 2 , IER ) 

CALL  CHECK ( I ER ) 

CALL  UNPACK ( 512 .TEMP, IN) 

DO  2  J-l ,20A8 

OUT ( J ) uCMPLX ( FLOAT (XN(J)) ,0,0) 

2  CONTINUE 

CALL  WRBLK( 2,32*1 , OUT , 32 , IER ) 

CALL  CHECK(XER) 

3  CONTINUE 

CALL  RESET 
STOP  VDTOCP 
END 


it  ’i’j’ji’jojiw.’r.i.v,’  av.v.’aa 
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c******* ••**•*•••*••***• *******  ******************************* 

c 

C  VCM.FR  -  DC  FORTRAN  5  -  LT  RICHARD  MILLS,  AUG  1984 
C 

C  Originally  CMULT  -  by  Cspt  Robert  Wells,  Dec  1983 

C  which  performed  coaplex  Multiplication  of  two  256  X 

C  256  complex  files  on  a  point  by  point  basis, 

c 

C  VCM  will  do  the  same  coaplax  Multiplication,  but  for 

C  two  coaplex  files  which  arc  N  x  N.  N  aust  be  a  power  of 

C  two  and  be  236  or  less  due  to  constraints  of  the  2DFFT 

C  prograas. 

C 

C  RELOAD  LINE: 

C  RLDR  VCM  0FLI1# 

C 

C  COMMAND  LIME: 

C  VCM 

C  Prograr  will  ask  for  input  and  output  file  naae,  create 

C  the  output  file  and  then  ask  for  file  size. 

C 

c«* ***************************** ****************************** 


COMPLEX  INI < 1024) , IN2 ( 1 024 ) , 0UT( 2 1 024 
INTEGER  INFLNM1 <  7  )  , INFLNM2 ( 7 ) , 0FLNM1 ( 7 ) 

ACCEPT  INPUT 

ACCEPT" ENTER  INPUT  FILENAME  #1  ->" 

READ(  1  1  ,  1  )  I.NFLNM1  (  1  ) 

ACCEPT" ENTER  INPUT  FILENAME  #2  ->" 

READ( 11,1) INFLNM2 ( 1 ) 

1  FORMAT ( S 1 3 ) 

ACCEPT" ENTER  OUTPUT  FILENAME  ->" 

READ( 11,1 )0FLNM1 ( 1 ) 

CALL  OPEMO  .  INFLNM1  ,  1  ,IER> 

CALL  CHECE(IER) 

CALL  OPEN( 2 , INFLNM2 , 1 , X ER ) 

CALL  CHECE(IER) 

CALL  DFILW  (0FLNM1.IER) 

CALL  CPILW  ( 0FLNM1 ,2 , KER ) 

CALL  CNECK(KER) 

OPEN  3,  0FLNM1 , ATT* "OR" , ERR* 1 00 

100  CONTINUE 
3  TYPE"  " 

TYPE"ENTER  SIZE  OF  FILES:" 

ACCEPT"  236,  128,  64,  OR  32  XSXZE 


* 


& 


<3 


C 

C  TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 
C 

I F( I S I ZE . EQ . 2 56 )  GO  TO  6 
IF(ISI2E.EQ.  128)  GO  TO  6 
IF(ISIZE. EQ .  64  )  GO  TO  6 
IF( ISIZE. EQ. 32 )  GO  TO  6 
GO  TO  5 

6  CONTINUE  ;  else  is  okay 

I FACTOR" 2 56/ISIZE  jeoapute  scaling  factor 

ITINES-64/ I FACTOR **2-1 

DO  3  I-O.ITIMES 

CALL  RDBLK( 1 . 16*1 , INI . 16, IER) 

CALL  CHECK(IER) 

CALL  RDBLK( 2t16*ItXN2t16tIER> 

CALL  CHECK(IER) 

DO  2  J" 1 , 102 A 

OUT(J)-INl(J)*CONJG(IN2(J)) 

2  CONTINUE 

CALL  VR>LK<3.16*I.0UT,l*.XBft> 

CALL  CHECK(IER) 

3  CONTINUE 

CALL  RESET 

STOP" <7X7X7  >VCH" 

END 


90 


o  o  a 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FVCM.FR  -  DC  FORTRAN  5  -  LT  RICHARD  MILLS,  19  SEPT  1984 

UPDATED  25  SEPT  1984 
UPDATED  3  OCT  1984 
UPDATED  4  OCT  1984 

FVCM  (  Fi  1 1  ered- Var  1  abl  e-Compl  ex-Mul  1 1  pi  1  cat  ion  )  is 
tlalltr  to  VCM  which  will  do  complex  mul tipi icat ion  of 
•ny  two  complex  files  which  are  H  x  N,  where  N  is  a 
power  of  two.  FVCM  differs  due  to  the  fact  that  it 
filters  out  everything  above  a  given  harmonic  out  of 
the  template  file.  This  program  anuses  that  the 
template  file  was  centered  at  the  origin  before  the 
2DFFT  was  taken. 

Presently  only  6A  x  84  and  236  x  256. 

RELOAD  LINE: 

RLDR  FVCM  9FLIB0 

COMMAND  LINE: 

TVCM 


COMPLEX  INI ( 1024) ,IN2< 1024) ,0UT( 1024) 
INTEGER  INFLNM1 (7) , INFLNM2 ( 7  > , OFLNMl ( 7 ) 


ACCEPT  INPUT 

ACCEPT"ENTER  INPUT  FILENAME  #1  ->" 
READ( 11,1) INFLNM 1(1)  •' 

ACCEPT "ENTER  I3PUT  FILENAME  #2  ->" 
READ( 1 1 , 1 ) INFLNM2 ( 1 ) 

1  FORMAT ( S 1 3 ) 

ACCEPT "ENTER  OUTPUT  FILENAME  ->" 
READ(11 , l)OFLNMl(l) 

CALL  0PEN( 1 , INFLNM1 , 1 , IER) 

CALL  CHECX(IER) 

CALL  0PEN( 2 , INFLNM 2 , 1 , I ER ) 

CALL  CRECK(IER) 

CALL  DFILV  (OFLNMl. IER) 

CALL  CFILW  (OFLNMl ,2, KER) 

CALL  CNECK(KER) 

OPEN  3,  OFLNMl , ATT-"0R" ,ERR» 1 00 
100  CONTINUE 
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non  non 


5  TYPE"  " 

TYPE"ENTER  SIZE  OF  FILES:" 
ACCEPT"  256,  or  64  — >",ISIZE 


C 

C  TEST  FOR  ONLY  THE  ALLOVARLE  SIZES 
C 

IF(ISIZE.EQ.25b)  CO  TO  7 
IF< ISIZE.EQ.64)  GO  TO  6 
GO  TO  5 

7  TYPE"  ”  ; dcttraln*  filter  (It*  for 

;  256  x  256  file. 

TYPE"ENTER  •i,-,32,  FOR  FILTER  OF  SIZE  7-235" 

TYPE"  (FILTER  SIZE  INCREMENTS  RY  S)  " 

TYPE"FOR  EXAMPLE:  i  FOR  A  7x7  FILTER,  2  FOR  A  15X15  FILTER," 
ACCEPT"  5  FOR  A  23X25  FILTER,  KTC...  — >  ",IFS 

IF( IFS . GT. 32 )CO  TO  7 
IFdFS.LT.  ;  )GO  TO  7 

IFSIZE*IFS*6- 1 
TYPE"  " 

TYPE"FILTER  TO  RE  USED  IS  ",IFSIZE,"  RY  " , IFSIZE , " . " 

TYPE"  * 

ACCEPT "ENTER  *  1 •  IF  TMIS  IS  CORRECT;  *0*  IF  MOT.  — >  ";IOX 
If < IOK -NE. 1 )G0  TO  7  ;rodo  input  if  incorrect 

6  CONTINUE  ;  else  ie  obey 

TYPE"-- - -  PROGRAM  RUNNING  - * 

I FACTOR*  2  56/ISIZE  ;co»pute  ecellnft  factor 

ITIMSS-64/IFACTOR**2-l 

FILL  OUTPUT  FILE  WITH  ZEROS 

DO  3  1*0 , ITIMES 

CO  2  J-1,1024 

OCT (J)*{0. 0,0.0)  ;fill  errey  OUT  with  teros 

2  CONTINUE 

CALL  WRRLR( 3,16*1, OUT ,16, IER )  ;write  errey  to  file 
CALL  CNKCK(IEt) 

3  CONTINUE  7 

IFCISIZE. EQ .4* )CO  TO  696  (file  ie  66x64 

MULTIPLY  AND  WRITE  TO  FILS  THE  DC  AMD  FIRST  THREE  HARMONICS 


u  u  u 


* 

4  *  a 


; do  first  block 


2  b 6  X  256 
699  I-32-IFS 


& 


CALL  RDBLK( 1  ,  16*1 , INI , 16, IER) 

CALL  CHECK (IER) 

CALL  RDBLK<2 , 16*1 , IN2 , 16, IER) 

CALL  CHECK ( I ER ) 

DO  18  J-l ,1024 

0UT( J )»( 0 . 0 , 0 . 0 )  ;  set  output  array  to  aaroi 

16  CONTINUE 

IVIDE-1 29-( ( IFSIZE-1 )/2 ) 

I START" 1 02 4- 3* I S I ZE+I WI DE 

{•tart  before  center  of  second  row 

XEND"ISTART+2*ISIZE 

;end  on  fourth  row 
DO  17  ICO" I START ,IEND,XSXZE 

{increment  by  ISXZE 

I STOP" IGO+I FS X  ZE- 1 

{do  XFSXZE  points  on  each| row 

DO  16  J"XG0 , XSTOP 

OUT( J)-XN1 (J)*CONJC(XN2(J)>  {Multiply  input  arrays 

16  CONTINUE 

17  CONTINUE 


.V 


s? 


A 

•V 

i 


i 


t 


CALL  URBLK (3,16*1, OUT ,16,1 ER) 
CALL  CHCCK(IER) 


'  •  **• 


C  COMPUTE  REMAINC  BLOCKS 
C 

IRBl-33-ir8 

IRB2-31+IFS 

DO  700  I"IRB1 , IRB2 


{do  reaalnlnR  blocks 


26 


CALL  RDBLK( 1,16*1, INI, 16,1 ER ) 
CALL  CHECK (IER) 

CALL  RDBLK(2, 16*1, IN2.16, IER) 
CALL  CHECK(IER) 

DO  28  J" 1 ,1024 

0UT(J)-<0. 0,0.0) 

CONTINUE 


I ST  ART" 1 VI DE 
I END- 1 START+3*  X  SIZE 


(•tart  before  center  of  first  row 
(end  on  fourth  row 
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i 


c, 

V 

«* 


ft 

ft 


I 

4 


V 

•j 


f 

r 


•WvV 

jaa jii  rf 


V  u  o 


00  27  ICO-ISTART,  HMD,  ISIZE 

IST0P-IG0+IFSIZE-1  ;do  IFSIZE  points  on  each  row 

DO  26  J "ICO , I  STOP 

0UT(J)*IN1(J)*C0NJG(IN2(J))  multiply  Input  arrays 

26  CONTINUE 

27  CONTINUE 

CALL  VRSLK(3, 16*1, OUT, 16.XER) 

CALL  CHECK(IER) 

700  CONTINUE 

GO  TO  900  ;skip  to  and 


6 A  X  6A  WITH  A  IS  X  IS  FILTER 


899  IMID3-(ITXMES+1 )/2  ;b!ock  past  cantor 

ZMID2-IMID3- 1  ; block  bafora  cantor 

TYPEMA  STANDARD  IS  X  IS  FILTER  IS  USED  ON  THE  64  X  64  FILE." 

I-IMID2  ; do  XMXD2  block 

CALL  RDRLXd.  16*1,  IN1.16.IER) 

CALL  CHECK(IER) 

CALL  RDRLR(2, 16*1 ,IN2, 16.IER) 

CALL  CHECKdES) 

DU  So 

OUT (J)*(0.0,0»0)  ;sat  output  array  to  saros 

58  CONTINUE 

I START* 1024-7* I S1ZE+27 

{Start  before  cantor  of  sacond  row 
{ 27*33{cantar )-< (lS-l)/2) 

I  END*  I START4-6*  1  SIZE  ;  and  on  sixtaanth  row 

DO  S7  ICO*ISTARi, isau.iSXZS  {increment  by  ISIZE 

I STOP* I 00+ 1 4  (do  IS  points  on  aach  row 

DO  S6  J*ICO, I STOP 

OUT  t  J ) • 1 N 1 ( J ) *C0NJC( I N2 ( J ) )  ;aultiply  input  arrays 

56  CONTINUE 

57  CONTINUE 

CALL  WR8LE( 5 ,16*1 ,  OUT, 16 , IER ) 

CALL  CHECKCER) 


•  4 


I-IMI03  ; do  IKID3  block 

CALL  RDBLK( I , 16*1 , INI , 16.IES 
CALL  CHECK  < 1 ER ) 

CALL  RDBLK( 2 , 1 6*1 , IN2 . 1 6 , IER) 

CALL  CHECK ( 1 ER ) 

DO  78  J»l  ,  1024 

OUT (J)*(0. 0,0.0/  ;i*t  output  array  to  itroi 

78  CONTINUE 

I START*2 7  ; itart  bafora  center  of  first  row 

IEND«ISTART-f7*ISIZE  ;«nd  on  eeventh  row 

DO  77  I GO* I START ,1END,ISXZE 

I  STOP*  I C04- 1  A  ;  do  ee*en  points  on  each  row 

DO  76  J*  I GO , I  STOP 

OUT ( J ) • I N 1 ( J ) *C0N JG ( X N2 ( J )  )  ;*ultiply  input  arrays 

76  CONTINUE 

77  CONTINUE 

CALI  VRILKO.  16*1  .OUT,  16, IER) 

CALL  CHECK(IER) 

m- 

»00  CALL  RESET 

STOP"  <  7X7X7  >FVCH" 

END 


A  A.V'  .V,1. 


•V/vV, 

■■ 1 


C  PROCRAM  CPTOVD  -  DC  FORTRAN  5  -  CAPT  ROBERT  WELLS,  DEC  1983 

C  THIS  PROGRAM  FINDS  1  HE  MAXIMUM  AND  MINIMUM  OF  A  256  X  256 
C  COMPLEX  FILF  AND  GIVES  THE  LOCATION  OF  THE  MAXIMUM.  ALSO,  THE 
C  PROGRAM  CONVERTS  THE  FILE  TO  PACKED  VIDEO  FORMAT  BY  PERFORMING 
C  A  LINEAR  SCALING  TO  THE  RANGE  OF  0  -  15  AND  TRUNCATING. 

INTEGER  IFLNM( 7 ) , OFLNM( 7 ) 

INTEGER  MA IN(7),F3(7),MS(2),S1(2),S2(2),S3(2) 

INTEGER  TEMP ( 5 12) ,OUT( 20*8) , I MAX , JMAX , MAX COL , MAXROV 
COMPLEX  IN ( 20*8 ) 

REAL  RMAG , MAXVAL .MINVAL, RANGE 

CALL  IOF( 2 , MA IN, IFLNM, OFLNM ,13, MS ,S1 , S 2 , S 3  ) 

CALL  OPEN( 1 , IFLHM, 1 , IER) 

CALL  CHECK  (I  Eft) 

CALL  OPEN! 2 . OFLNM , 3 , I ER ) 

CALL  CHECK ( I ER ) 

MAXVAL-0.0 

MINVAL-1E50 

DO  3  1-0,31  ;  FIND  MAX  AND  MIN 

CALL  RDBLKf 1 ,32*1, IN, 32, I ER ) 

CALL  CHECK(IER) 

DO  3  J- 1 , 20*8 

RMAG -CABS ( I N ( J ) ) 

I F ( RMAC . LE . MAXVAL) GO  TO  2 
MAX VAL- RMAG 


IMAX-I 
JMAX— J 

2  I F<  RMAC . IT.MINVAL)M 

3  CONTINUE 
RANGE-MAX VAL-M1NVAL 


snvalI 


RMAC 


MAXCOL-NOD< JMAX-! ,256)* 1 
MAXROV^*  *  IMAX<MMT<<  JMAX-1  )/!«34)M 


TYPE-MAX 

COL 

• 

- , MAXCOL 

TYPE-MAX 

ROW 

• 

”, MAXROV 

TYPE-MAX 

VAL 

m 

- .MAXVAL 

TYPE-MIN 

VAL 

m 

*, MINVAL 

DO  5  I«0,3i 


CALL  R03LK{ 1 ,32*1 .IN , 32 , IER) 
CALL  CHECKOER) 


A  J«1 , 2046 

Jr!.<,i>“AN1NT(15*0*(CABS(IN(J>)-MlNVAL)/RANGE) 


CONTINUE 


CALL  REPACK(5J2 .OUT, TEMP) 
CALL  WRBLK'2.2*I,TEMP,2,IER) 
CALI  CHECK(IER) 

5  CONTINUE 


CALL  RESET  f.  , 

STOP0<7><7><7X7>CPtOVDm 

END  ’ 


non 


c 

C  SMALL.FR  -DC  FORTRAN  3-  LT.  RICHARD  L.  MILLS  6  AUG  1964 
C 

C  This  program  will  accept  the  name  of  a  video  file  and 
C  create  a  new,  smaller,  file  which  contains  only  specific 

C  video  info.  The  only  other  info  needed  is  the  X-Y 

C  coordinates  (Row-Column)  of  the  original  position  of 
C  interest  and  the  size  of  the  area  to  be  copied  into  the 
C  smaller  file. 

C 

C  RELOAD  LINE: 

C  RLDR  SMALL  F3PICBUF.LB  0FLIR0 

C 

C  COMMAND  LINE: 

C  SMALL 

C 

c* ***************************** ft*************** **************** 

C*  PICBUF  STUFF 

PARAMETER  NBUFSZ-300 

PARAMETER  NCOLSMAX*236 

INTEGER  I BUF ( NBUFSZ ) , IBUF2 (NBUFSZ) 

INTEGER  IARRAY ( NCOLSMAX ) , INAME(20) , IOTNM( 20) 

IMPLICIT  TNTEGER (A,B,C,R,X,Y) 

C*  I/O  CONSTANTS 

PARAMETER  IN- 11 
PARAMETER  OUT- 10 
C 

C  ACCEPT  INPUTS 
C 

ACCEPT” NAME  OF  THE  INPUT  FILE: — >" 

READ ( IN , AO ) INANE ( 1 ) 

40  FORMAT ( S40 ) 

ACCEPT" NAME  OF  THE  OUTPUT  FILE: — >" 

READ( IN, 41 ) I0TNM( 1 ) 

41  FORMAT(SAO) 

GET  SIZE  OF  NEW  IMAGE 


*  *  *  «• 

TYPE"ENTER  SIZE  OF  AREA  OF  INTEREST:” 
ACCEPT"  12®,  64,  32:  ”,ISIZE 
C 

C  TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 

C  (DUE  TO  CON5TRANTS  ON  FFT  A  I FFT  PROGRAMS) 

C 

IF(ISIZE.EQ. 128)  GO  TO  6 
IF(ISIZE. EQ .64)  CO  TO  6 
IF(ISIZE.EQ. 32)  GO  TO  6 
GO  TO  5 

CONTINUE  ;  size  is  okay 


6 


o  o  o 


c 

c 

C 


C 

c 

c 


I  FACTOR- 2 5 6 / I S I ZE  {calculate  reduction  factor 

GET  POSITION 
Type"  " 

TYPE"ROW, COLUMN  OF  THE  POINT  OF  INTERST  * 

' ACCEPT" IN  THE  ORIGINAL  IMAGE — >*,R,C 
X-C 
Y-R 

PICBUFF  FOR  INPUT  FILE 

CALL  P1CFMT( IBUF , 1NAME , IK SR) 

CALL  CFMT(  IBUF  ,NROWS  ,t!COLS, MBITS, XMODE) 

CALL  MAKB(IBUF) 

CaLL  PICIN(IBUF, INAME, IHDR) 


CALL  DFILW  ( IOTNM, XER)  {delete  output  file  if  it 

IF< IER.NE. 1 )GOTO  90 
90  CONTINUE 

PICBUFF  FOR  OUTPUT  FILE 


CALL  GFMT( I BUF , NROWS , NCOLS , NBITS , IMODE ) 
NROWS2-NROWS/II ACTOR 
NCOLS 2* NCOLS/ 1 FACTOR 

CALL  PFMTORUF2  , NROVS2 , NC0LS2 , NBITS . IMODE > 
CALL  MAKB ( I BUF2 ) 

C 

C  A INSURE  PROPER  DISTANCE  FROM  BORDERS 
C 

IDIST-ISIZE/2 
IF(IDIST.LE.R)CO  TO  20 
R-IDIST 
TYPE"  " 

TYPE"  TOO  CLOSE  TO  IMAGE  BORDER," 

WRITE  FREE(OUT)"  NEW  ROW  IS  ",R 

20  ,IF(IDIST.LE.C)GO  TO  21 
C-IDIST 

TYPE"  " 

TYPE"  TOO  CLOSE  TO  IMAGE  BORDER," 

WRITE  FREE (OUT) "NEW  COLUMN  IS  ",C 

21  IF( ( 256-R) .GE.ZDXST)GO  TO  22 
R»256-IDIST 

TYPE"  " 

TYPE"  TOO  CLOSE  TO  IMAGE  BORDER," 

WRITE  FREE(OUT)"  NEW  ROW  IS  ",R 

22  IF( ( 2S6-C) . GE. XDXST)GO  TO  23 
C-2SS-IDXST 


exiet 
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TYPE"  " 

TYPE"  TOO  CLOSE  TO  IMAGE  BORDER," 

WRITE  FREE  ( OUT ) "NEW  COLUMN  IS  ",C 
23  CONTINUE  ;R,C  arc  at  least  IDIST  fros  isagc  border 

TYPE"  " 

TV  PE"MOV I NG  PIXELS" 

1START-R-1DIST 
I  STOP" I START+1 S IZE- 1 . 

IROW2-1 

ICOL-C-IDIST 

ICOL2-1 

DO  200  I ROW" I  START . X  STOP 

CALL  GROW< IIUF, IROW, ICOL, ISIZE, I ARRAY) 

CALL  PROW( I BUF2 , IROW2 , XCOL2 , I SIZE , I ARRAY ) 

I ROW 2 -IROW2+1 
200  CONTINUE 


CALL  PICOUT( IBUF2 , lOTNM, IHDR)  (write  to  file 

CALL  RELB(IBUF)  ;close  buffers 

CALL  RELB( IBUF2 ) 

CALL  RESET 

STOP" <7 ><7><7><7> SMALL" 

END 


100 


rtr.fc-  wr.- 


non 


c* ******************************************************** ****<< 

c 

C  VSPIN.FR  -DC  FORTRAN  5-  LT  RICHARD  L.  MILLS  17  AUG  1984 
C  UPDATED  7  SEPT  1984 

C 

C  Thi*  program  will  accept  the  name  of  a  file  which  contains 
C  a  template  and  create  a  new  file  which  places  the  template 

C  center  at  the  origin  of  the  video  file.  This  is  much  like 

C  the  program  SPIN.FR  except  VSPIN.FR  ( vari able-S pin > works 

C  on  variable  sized  video  files  and  creates  files  of  the 
C  same  size.  The  input  is  assumed  to  bo  cantered,  there- 
C  fore  no  coordinates  are  needed. 

C 

C  RELOAD  LINE: 

C  RLDR  VSPIN  FSPICBUF.LB  0FLIB0 

C 

C  COMMAND  LINE: 

C  VSPIN 

C 

c*** *******************  **************************************** 

C*  PICBUF  STUFF 

PARAMETER  NBUFS2-300 

PARAMETER  NCOLSMAX-256 

INTEGER  IBUF(NBUFSZ) , I BUF2 (NBUFSZ) 

INTEGER  1 ARRAY ( NCOLSMAX ) , I NAME ( 20) , IOTNM( 20 ) 

IMPLICIT  INTEGER (A,B,C,R,X,Y) 

C*  I/O  CONSTANTS 

PARAMETER  IN- 11 
PARAMETER  OUT- 1 0 
C 

C  ACCEPT  INPUT 
C 

ACCEPT  "NAME  OF  THE  INPUT  FILE:— >" 

READ(IN,40) I NAME ( 1 ) 

40  FORMAT(S40) 

ACCEPT "NAME  OF  THE  OUTPUT  FILE: — >" 

READ( IN , 4 1 ) I0TNM( 1 ) 

41  FORMAT ( S40 ) 

5  TYPE"  " 

TYPE"ENTER  SIZE  OF  INPUT  FILE:" 

ACCEPT"  256,  128,  64,  OR  32  — >".ISIZE 

TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 

IF< ISIZE.EQ. 256)  CO  TO  6 
IF(ISIZE.EQ. 128)  CO  TO  6 
IF(ISIZE.EQ.64)  GO  TO  6 
IF( ISIZE.EQ. 32)  CO  TO  6 
GO  TO  5 

6  CONTINUE  ;  sis*  is  okay 
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c 

C  PICBUFF  FOR  INPUT  FILE 
C 

NROVS-ISIZE 

NCOLS-XSIZE 

NBITS-4 

INODE-1 

IHDR-1 

CALL  PFMT( IBUF.NROUS, NCOLS, NBITS, INODE) 

CALL  NAKB(IBUF) 

CALL  PXC1N(XBUF, I NAME  . XHDR ) 

CALL  DFILW  (IOTNM.IER)  ;d«l«t«  output  ill*  If  it 

IF( IER.NE. 1 )GOTO  90 
90  CONTINUE 
C 

C  PICBUFF  FOR  OUTPUT  FILE 
C 

CALL  GFMTC 1BUF.NROVS , NCOLS , NBITS . IKODE ) 

CALL  PFMT( I BUF2 , NROWS , NCOLS , NBITS , 1 MODE ) 

CALL  MAKB(IBUF2) 

C 

C  CALCULATE  MID  PT. ,  QUARTERED  SIZE,  CENTER  STARTING  PT. 

C 

XM-I SIZE/2 
YM-ISIZE/2 
ILEN-ISI ZE/2 
I  START* ( IS IZE/ 2 )+ 1 
C 

C  MOVE  UPPER  LEFT  CORNER 
C 

TYPE"MOVING  PIXELS” 

IROW2- YM+1 

ICOL-1 

ICOL2-XM+1 

DO  200  IROW* 1,1 LEN 

CALL  CROW(  XBUF  1 1  ROW’,  ICOL.  I  LEN  . 1  ARRAY) 

CALL  PROW( XBUF2 , XROV2  , ICOL2 . ILEN, IARRAY) 
IROW2-IROW2+ 1 
200  CONTINUE 
C 

C  MOVE  UPPER  RIGHT  CORNER 
C 

TYPE"MGVINC  PIXELS" 

1ROV2-YM+I 
XCOL-XM+1 
ICOL2- 1 

DO  400  IROW- 1,1 LEN 

CALL  CROV( IBUF, IROW, ICOL, ILEN, IARRAY) 

CALL  PROV(IBUF2,XROV2,XCOL2. ILEN, IARRAY) 
IROV2-IROW2+1 
400  CONTINUE 


102 


C  MOVE  LOVER  LEFT  CORNER 
C 

TYPE"MOVING  PIXELS" 

XROW2-1 
I COL" I 
XCOL2-XM+1 

00  300  I  ROW* I  START, ISIZE 

CALL  GROW( IRUF, I ROW , I COL . I LEN , I ARRAY ) 

CALL  PROW ( I SUF2 , XR0W2 , XC0L2 , I LEN , I ARRAY ) 
IROW2-IROV2+1 
300  CONTINUE 
C 

C  MOVE  LOWER  RIGHT  CORNER 
C 

TYPE"NOVING  PIXELS" 

IR0V2-1 

ICOL-YM+1 

IC0L2-1 

DO  500  IROV«ISTART, IS1ZE 

CALL  CROV( IRUF , I ROW, I COL, I LEN 4  X ARRAY) 

CALL  PROV( 1RUF2 , IR0W2 , IC0L2 , XLEN , I ARRAY ) 
IROV2-IROV2+1 
500  CONTINUE 

CALL  PICOUT( IRUF 2 » IOTNM, 1 )  ;wrlt«  to  file 

CALL  RELR(IRUF)  jcIom  buffers 

CALL  RELB( ISUF2 ) 

CALL  RESET 

8TOP"<7><7X7><7>VSPIH" 

END 
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c* ******* ft* ***************************************************** 

C 

C  NORM . FR  -  DC  FORTRAN  5  -  LT  RICHARD  MILLS,  11  SEPT  1984 

C 

C  This  prograa  will  accept  a  template  file  and  a  target 

C  file  and  create  a  third  file  which  will  be  a  new  target 

C  file  with  it*  energy  normalized  to  that  of  the  template. 

C  File  sizes  may  vary;  32,  64,  128,  or  256.  (Normalized 

C  energy  may  not  exactly  equal  that  of  the  template  due 

C  to  round  off  errors.) 

C  This  program  is  similar  to  N0RM2.FR  except  for  the  fact 

C  NORM . FR  will  automatically  continue  the  normalization 

C  process  repetitively  until  it  can  no  longer  do  ao  due 

C  to  the  round  off  error. 

C 

C  RELOAD  LINE: 

C  RLDR  NORM  UNPACK  REPACK  0FLI80 

C 

C  COMMAND  LINE: 

C  NORM 

C 

c***«* **********************a*********************************** 


INTEGER  TEMP(2S6) , I N ( 1024) 

INTEGER  I NFLNM( 7 ) , 0FLNM1 ( 7 ) , 0FLNM2 ( 7 ) 
REAL  OUT <  1024) 

DOUBLE  PRECISION  ENERGY  1  .F.NERGY2  .ENERCY3 
ACCEPT  INPUT 


ACCEPT"ENTER  TEMPLATE  FILENAME  ->" 

READ< l 1 , 1 )INFLNM< 1 ) 

1  FORMAT ( S 1 3 ) 

ACCEPT" ENTER  TARGET  FILENAME  ->" 

READ ( 11,1 )OFLNMl ( 1 ) 

ACCEPT "ENTER  NEW  NAME  OF  NORMALIZED  TARGET  •>" 
READC1 1 , 1 )0FLNM2( 1 ) 

CALL  OPEN ( 1 , INFLNM , 1 , IER ) 

CALL  CHECK(IER) 

CALL  OPEN (  2 , 0FLNM.1  ,  1  , 1 ER  ) 

CALL  CHECK ( IER) 

CALL  DFILW  ( OFLNM2 , I ER ) 

CALL  CFILW  ( 0FLNM2 , 2 , XER  ) 

CALL  CHECK(KER) 

OPEN  3,  0FLNM2 , ATT* "OR" , ERR* 100 
100  CONTINUE 
5  TYPE"  " 

TYPE "ENTER  SIZE  OF  FILES:" 

ACCEPT"  256,  128,  64,  OR  32  — >"iISIZE 
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c 

C  TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 
C 

IF< ISIZE.EQ. 256)  GO  TO  6 
IF( ISIZE.EQ. 1 26)  GO  TO  6 
IF(ISIZE.EQ.6A)  GO  TO  6 
IF(ISIZE.EQ.32)  CO  TO  6 
GO  TO  5 

6  CONTINUE  ;  six*  la  okay 

I  FACTOR-2  5 6 / 1 S I ZE  {calculate  reduction  factor 

IT I MBS *64/ I FACTOR** 2-  1 

ENERGY  1  *0  {Initialise  energies  to  itero 

ENERGY2-0 
C 

C  REAO  TEMPLATE  ENERGY 
C 

DO  3  I-O.ITI'JZS 

CALL  RDBLK ( 1 , I , TEMP , 1 , IER) 

CALL  CHECK(IER) 

CALL  UN PACK ( 236 .TEMP , IN) 

DO  2  J- 1.1 02* 

ENERGY  1  •ENERGYl+t  IN(J))**2 

2  CONTINUE 

3  CONTINUE 
C 

C  READ  TARGET  ENERGY 

C 

DO  13  I-O.ITIMES 

CALL  RDRLK(2 , I .TEMP. 1 , IER) 

CALL  CHECK(IER) 

CALL  UN PACK ( 236. TEMP. IN) 

DO  12  J-1.10?a 

ENERCY  2 "ENERGY 2  +  (  IN(J))*"2 

12  CONTINUE 

13  CONTINUE 
ENERGY3*ENERGY2 

DISPLAY  DATA 

TYPE"  " 

TYPE "TEMPLATE  ENERGY  j  ". ENERGY! 

TYPE"ORICINAL  TARGET  ENERGY  :  " , ENERCY2 
30  RATI0*DSQRT(ENERGY1/ENCRGY2) 

TYPE"S0.  RT.  ENERGY  RATIO  :  ".RATIO 
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READ  TARGET  ENERGY  AND  NORMALIZE 

ENERGY2-0 
DO  23  I  »0  , ITI MES 

CALL  RDBLK(2.I .TEMP, 1 ,IER> 

CALL  CHECK ( I ER ) 

CALL  UKPACK(236 , TEMP , IN) 

DO  22  J-1,1024 

OUT< J)-IN( J) •RATIO 

IN ( J ) • I P IX ( OUT ( J )+0 . 5 )  {round  off 
IF(IN(J).CT.1S)  IN(J)-15 
ENERGY 2 "ENERGY 2+( XN(J))**2 

22  CONTINUE 
CALL  REPACK ( 236, IN, TEMP) 

CALL  VRRLK(3,I  .TEMP,  I  ,IER) 

CALL  CHECK(XER) 

23  CONTINUE 

TYPEMNEW  TARGET  ENERGY  :  " a  ENERGY 2 

TEST  IF  ENERGY  IS  AS  CLOSE  AS  IT  CAN  GET 

IF (ENERGY 3 . EQ . ENERGY 2 )C0  TO  31 
ENERGY 3* ENERGY 2 

CALL  CLOSE(2,IER) 

CALL  CHECK(XER) 

CALL  DFILV  ( "HOLD . VD* , I ER ) 

CALL  CFILW  ("HOLD. VD",2,KER) 

CALL  CHECK(KER) 

OPEN  A,  "HOLD.VD",ATT»"OR*»ERR»!01 
1  CONTINUE 

TRANSFER  OUTPUT  TO  HOLD. VD  AND  CHANGE  CHANNELS 

DO  33  I *0 , ITXMES 

CALL  RDtLK(3,I .TEMP, 1 .IER) 

CALL  CHECK ( IER ) 

CALL  VRSLK(A,X, TEMP, 1,IER) 

CALL  CHECK(XER) 

33  CONTINUE 

CALL  CL0SE( A , IER ) 

CALL  CHECK(XER) 

CALL  OPEN (  2  , "HOLD . VD", 1 . IER ) 

CALL  CHECK ( IER ) 

GO  TO  30 
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5 1  CALL  RESET 

CALL  DFILW("HOLD.VD",IER) 
CALL  CHECK ( I ER  ) 
STOP"<7><7><7><7>NOKN" 

End 


N* 

£ 

L 


I*' 


t 

C: 

V/ 

r; 

r. 
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C *  ft**************************  *******  ***************  ******  ****** 

c 

C  VTC.FR  -  pc  FORTRAN  3  -  LT  RICKARD  HILLS,  AUG  1984 
C 

C  Originally  VDTOCP  -  by  Capt  Robart  Walla,  Dac  1983 

C  which  converted  a  236  X  236  packed  vadlo  file  to  a 

C  236  X  236  coaplex  file. 

C 

C  VTC  will  de  the  saae  coaplex  file  conctruction  but  for 

C  any  packed  video  file  which  ie  N  x  N.  N  suit  be  a  power 

C  of  two  and  be  236  or  leee  due  to  conetrante  on  the  2DFFT 

C  prograae. 

C 

C  RELOAD  LINE: 

C  RLDR  VTC  UNPACK  fFLIlf 

C 

C  COMMAND  LINE: 

C  VTC 

C  Prograa  will  aak  for  input  and  output  file  naae,  create 

C  the  output  file  and  then  aak  for  file  alee. 

C 

Ceeeeeeeeeeeeeaeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeaeei 

INTEGER  TEMP( 236),IN(1024) 

INTEGER  INFLNM(  7 ) , OFLNMl (  7  ) 

COMPLEX  OUT( 2048 ) 

C 

C  GET  INPUTS 
C 

ACCEPT"ENTER  INPUT  FILENAME  ->" 

READ( 1 1 , 1 > I NPLNN( 1 ) 

1  FORMAT ( S 1 3 ) 

ACCEPT "ENTER  OUTPUT  FILENAME  ->" 

READ( 11,1 )0FLNM1 ( 1 ) 

CALL  0PEN( 1 , INFLNN, 1 , IBS) 

CALL  CHECX(IER) 

CALL  DFILV  (OFLNMl, IBS) 

CALL  CFILV  (OFLNMl ,2, KER) 

CALL  CHECK(KER) 

OPEN  2,  OFLNMl ,ATT-"OR",ERR»lOO 
100  CONTINUE 
C 

C  GET  FILE  SIZE 

C 

3  TTPE"  " 

TYPE"ENTER  SIZE  OF  INPUT  FILE:* 

ACCEPT"  236,  12S,  64,  OR  32  — >".ISIZE 
C 

C  TEST  FOR  ONLT  THE  ALLOWABLE  SIZES 
C 


6 


C 

C 

c 


IF< ISIZE.EQ.256)  CO  TO  6 
IF( TSIZE. EQ. 1 28)  CO  TO  6 
IFOSIZE.EQ.64)  CO  TO  6 
IFOSI2E.EQ.32)  GO  TO  6 
GO  TO  5 

CONTINUE  ;  alse  la  okay 

I  FACTOR-2  56/ISIZE 
ITIMES«64/ I FACTOR* *2- I 

DO  CONVERSION  TO  COMPLEX 

DO  3  I-O.ITIMES 


;|tt  scaling  factor 


CALL  RDBLK( 1 , I .TEMP , 1 , IER) 
CALL  CNECK(IER) 

CALL  UN PACK ( 236, TEMP, ZN) 


2 


DO  2  J-1,1024 

OUT ( J ) -CMP LX (FLOAT  <IN(J))tO.O) 

(act  imaginary  part  to  zaro 

CONTINUE. 


CALL  VRBLK(2,14*I,OUT.16.ZER) 
CALL  CHECK ( I ER ) 


3  CONTINUE 
CALL  RESET 

STOP"<7><7><7  ><7> VTC" 

END 


y 

A* 


w 
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c************************************************************* 

c 

C  CTV.FR  -  DC  FORTRAN  5  -  LT  RICHARD  MILLS,  AUC  1984 
C 

C  Originally  CPTOVD  -  by  Capt  Robert  Welle,  Dec  1983 

C  which  converted  a  256  X  256  coaplex  file  to  a  256  X  256 

C  packed  video  file. 

C 

C  CTV  will  do  the  i«et  coaplex  file  conetructloh  but  for 

C  and  packed  video  file  which  ie  N  x  N.  N  auat  be  a  power 

C  of  two  and  be  256  or  lose  due  to  conatrante  on  the  2DFFT 

C  prograas. 

C 

C  RELOAD  LINE: 

C  RLDRCTV  IOF  REPACK  «FLIS9 

C 

C  COMMAND  LINE: 

C  CTV 

C  Prograa  will  aek  for  input  and  output  file  nano,  create 

C  the  output  file  and  then  aek  for  file  (lit, 

C 

C*ee*eeeee*6eeeee*eeaeeeeeaaeaaeece*e*eaeeeaaeeeeeeeeeeeeeeeee< 
INTECER  I FLNM( 7 ) , OFLNM( 7 ) 

INTEGER  TEMP( 2  56 ) , OUT( 1024), IMAX , JMAX , MAXCOL , NARROW 
COMPLEX  iN(  1024) 

REAL  RMAG , MAX VAL . MI NVAL , RANGE 

ACCEPT  INPUT 

ACCEPT"ENTER  INPUT  FILENAME  ->" 

READ( 1 1 , 1 ) 1FLNM( 1 ) 

1  FORMAT ( S 1 3 ) 

ACCEPT1* ENTER  OUTPUT  FILENAME  ->" 

READ( 11,1 )OFLNM( 1 ) 

CALL  OPEN( 1 , IFLNM, 1 , IER) 

CALL  CHSCK(IER) 


CALL  DFILW  (OFLNN, IER) 

CALL  CFILW  (OFLNM, 2 , XER) 

CALL  CHECK(KER) 

OPEN  2,  OFLNM, ATT- "OR", ERR- 100 


CONTINUE 
TYPE"  * 

TTPE-INTRR  SIZE  OF  INPUT  FILE:" 

ACCEPT"  256,  128,  64,  OR  32  — >",ISIZE 

TEST  FOR  ONLY  THE  ALLOVASLE  SIZES 


6 


C 


IF( I  SIZE . EQ .  2  56  )  GO  TO  6 
IF( ISIZE.EQ. 128)  GO  TO  6 
IF(  ISIZE . EQ. 64  )  CO  TO  6 
1F( ISIZE.EQ. 32)  GO  TO  6 
CO  TO  5 

CONTINUE  i  llte  Is  okay 

IFACTOR-256/IS I2E 
I TIMES- 64/ I  F  ACTOR* *2-1 

MAX VAL-0 . 0 
MINVAL-1E50 


;coaputc  scaling  factor 


;aat  initial  valuci 


00  3  1-0 , ITINES  ; f ind  aax  and  ain 

CALL  RDSLK( 1 , 16*1 , IN, 16, IER) 

CALL  CHECK ( I  EE ) 

DO  3  J-1,1024 

■  MAG-CABS ( I N ( J  )  ) 

I F ( KMAG .LE.MAXVAL)GO  TO  2 
MAX V AL-RMAG 
IMAX-I 
JMAX- J 

2  I F( RMAC . LT .MI MVAL )M I NVAL-BMAG 

3  CONTINUE 


RANGE-MAX VAL-M I N VAL  ;find  rang*  of  valuaa 

MAX COL- MO D ( JMAX-1 ,  2  5  6  )  +  1 
MAXROV-B  *  I  MAX-*’  I  NT  (  (  JNAX-  1  )/256)-fl 

TYFE"MAX  COL  -  ".MAXCOL  ;  REFORT  MAX  INFO. 

TTFE"MAX  ROW  -  ".MAXROW 

TYFE*MAX  VAL  -  " , MAXVAL 
TYFE"MIN  VAL  -  * , MINVAL 

DO  I  1-0 , 1  TIMES  ;r«write  aa  vidao 

CALL  RDBLK( 1 , 16*1 , IK. I  6, IER) 

CALL  CHECK(XER) 

DO  4  J-l . 1024 

OUT< J)-ANINT( l  5.0* <CABS<IN<J ))-MINVAL)/RANCE) 

4  CONTINUE 


1  1  1 


CALL  REPACK( 256 , OUT, TEMP) 
CALL  VRBLK(2 , I .TEMP, 1 , IER) 
CALL  CHECK(IER) 

S  CONTINUE 

CALL  RESET 

STOP"<7X7X7X7>CTV" 


£*•***********••••••**»****•*•***********•***•••**•••*••**••*• 

c 

C  EVAL.FR  -  DC  FORTRAN  5  -  LT  RICHARD  MILLS,  12  SEPT  1944 

C  Updated  17  SEPT  1984 

C 

C  This  program  will  evaluate  the  correlation  peaks  after 

C  normal i cat  ion .  The  evaluation  is  based  on: 

C 

C  1)  Peak  location 

C  2)  Peak  shar pnes s / uniqueness . 

C 

C  A  hard  copy  output  may  be  obtained  by  printing  the 

C  file  “LPT" . 

C 

C 

C  RELOAD  LINE: 

C  RLDR  EVAL  F5PICBUF.LB  0FLXB0 

C 

C  COMMAND  LINE: 

C  EVAL 

C 

c* A*********************************************************** 

C  PICBUF  STUFF 

PARAMETER  NBUFS2-300 
PARAMETER  NCOLSMAX-S12 
INTEGER  I BUFF3 ( NBUFS2 ) 

INTEGER  I  ARRAY ( NCOLSMAX ) 

C* 

PARAMETER  NLIMIT-20 
PARAMETER  PERCENT*  5.0 
INTEGER  TEMP (236},IN(1024) 

INTEGER  MA I NFL ( 7 ) 

INTEGER  INFLNM(400) 

INTEGER  ISI2E(NLIMIT) 

REAL  DIST(NLXMIT) 

REAL  MIDPT 

INTEGER  MAXROV(NLIMIT) , MAXCOL< NLXMXT) 

REAL  SCORE  l  (NLXMIT)  .SC0RF.2  (NLXMZT)  ,  IRANK(NLXMXT) 

REAL  HIS( NLXMIT, 16) 

CALL  DFILW  (”LPT”,XER)  ;open  up  output  file 

OPEN  12, "LPT" 

C 

C  SET  "HIS"  ARRAY  TO  2ER0S 
C 

DO'  40  1-1  ,  16 

DO  41  ICOUNT-1 .NLIMIT 
H I S ( I COUNT , I ) *0 . 0 
41  CONTINUE 

40  CONTINUE 
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ICOUNT-O 


set  counter 


ACCEPT  INPUTS 
Type"  ■ 

ACCEPT" ENTER  FILENAME  OF  THE  MAIN  FILE  ->" 

READ ( 1 1 , 1 )MAINFL( 1 ) 

I  COUNT*  I  COUNT*’  1  {Increment  counter 

TYPE"  " 

ACCEPT "ENTER  FILENAME  ->" 

READ (  1  1  ,  I  )INFLNN(  ICOUNTMO) 

1  FORMAT( S20) 

CONTINUE 

ACCEPT" ENTER  SIZE  OF  FILE:  256,  128,  6A,  OR  32  — >", 
♦  ISIZE(  I  COUNT) 

TEST  FOR  ONLY  THE  ALLOWABLE  SIZES  DUE  TO  FFT  PROGRAMS 

IF(ISIZE( I COUNT ) . EQ .256)  GO  TO  6 
IF(ISIZE(ICOUNT) .EQ. 128)  GO  TO  6 
IF( I S I ZE( 1 COUNT ) . EQ .64)  GO  TO  6 
IF( ISIZE( I COUNT ).EQ.32)  60  TO  6 
GO  TO  5 

CONTINUE  ;  else  is  okey 

IF ( I COUNT . EQ . N LI NIT) GO  TO  8 
TYPE  ICOUNT,”  FLIES  SO  FAR" 

ACCEPT  "ENTER  1  FOR  ANOTHER;  0  TO  EVALUATE  ".IMORE 
IF( IMORE.EQ. 1 )C0  TO  7 
TYPE"  " 

TYPE" EVALUATING" 

TYPE"  " 

NUMFILES- ICOUNT 
INPUT  COMPLETE,  NOV  EVALUATION 


DO  9  ICOUNT- 1 .NUMFILES 

CALL  OPEN ( 1 , INFLNN( ICOUNT* 10) , 1 ,IER) 
CALL  CHECK(IER) 

CALCULATE. MID  PT. 

MI DPT- I SIZE (I COUNT)/ 2+0.5 

PICBUFF  FOR  INPUT  FILE 

NROWS-I SIZE 
NCOLS-ISIZE 
NBITS-4 
INODE- 1 


1 1  A 


IHDR- 1 

CALL  PFMT ( I BUF , NROWS ,  NCOLS.NBITS .INODE) 
CALL  HAKB(IBUF) 

CALL  PICIN( I BUF. INFLNM( 1  COUNT*! 0) .IHDR) 


C  FIND  MAX  LOCATION 
C 

IFLAG-0 

MAXROV( ICOUNT)-999  (inltallly  set  to  vary  high  values 

MAXCOL( I COUNT ) *999 
DO  20  IROW-1 ,  I S 1 ZE ( I  COUNT ) 

DO  21  ICOL- 1 , I S I ZE ( ICOUNT) 

CALL  GPNT ( I BUF , I  ROW , I COL ,  X VAL) 

MIS (ICOUNT, I VAL+1)-HIS( ICOUNT, XVAL+1 )♦!.  0 

(construct  histrograa  data  array 
IF( I FLAG . EQ . 1 )GO  TO  21  (skip  If  found 

IF(I VAL.EQ. 1 5 > IFLAG-1  ;flnd  sax  point 

I  F(  I  FLAG  .  EQ  .  1  )MAXROV(  ICOUNT)-IROV 
I F ( 1  FLAG . EQ . 1 ) MAXCOL ( I COUNT )- ICOL 
21  CONTINUE 

20  CONTINUE 

CALL  RELB(IBUF) 


i  ® 


C  EVALUATE  MAX  LOCATION 
C 

01 ST( I COUNT) -SQRT( (MI DPT-MAXCOL( ICOUNT) )**2+ 

.  ♦  (MIDPT-MAXROV( ICOUNT) )**2) 

SCORE  1 (ICOUNT)- 1 .0-(2*DIST( XCOUNT)/ISIZE( ICOUNT)) 

I F ( SCORE 1 ( ICOUNT) . LT .0.0)5 C0REI( ICOUNT) -0 . 0 
C 

C  DETERMINE  SHARPNESS /UNIQUENESS  OF  PEAK 
C 

PERCENTAREA- ( PERCENT/ 1 00 . 0 ) * I S I  ZE ( I COUNT ) *  *  Z 
I SUM  1-0 
ISUN2-0 
J- 1  6 

25  ISUH1-ISUM1+KIS( ICOUNT, J) 

I F( I SUM 1 . GT . PERCENTAREA ) GO  TO  26  ("paakness"  dataralnad 

ISUM2-ISUM1 

J-J-l 

GO  TO  25 

26  FRACTION- I SUM2 /PERCENTAREA 

S POT- J+ FRACTION  (deternlne  score!  function's  input 


C  EVALUATE  THE  SHARPNESS/UNIQUENESS 
C 

SCORE 2 ( ICOUNT) -SQRT ( ( 1 6 . 0-S POT ) / 1 0 . 0 ) 
IF(SCORE2( ICOUNT). GT. 1 . 0 ) SC0RE2 ( ICOUNT) • 1 . 0 

SCORE ( ICOUNT) -SCORE  1 (ICOUNT) *SC0RE2( ICOUNT) 


u  o  o 


I RANK ( I  COUNT) -I FIX( SCORE ( I COUNT) M 0 )♦ 1 
IF( I RANK ( I  COUNT) . GT . 1 0 ) I RANK ( I COUNT) » 1 0 

CALL  CLOSE ( 1 . IER) 

CALL  CHECK ( I ER ) 

9  CONTINUE  ; END  OF  EVALUATION  LOOP 

WRITE  (10, lA)MAINFL(l) 

WRITE  FREE  (10)  "FILE  NAME  SCORE1  SCORE2 

♦  TOTAL  SCORE  RANK" 

WRITE  FREE  (10)  " -  -  - 


DO  AS  ICOUNT-1 ,NUMFILES 
ICNT»ICOUNT  ;  to  uv«  coo* 

I TEN «I COUNT* 1 0 

WR ITE( 10,2) INFLNM( I TEN ) , SCORE 1 ( 2 CRT) , SCORE 2 ( I CNT ) , 

♦SCORE ( ICNT) , I RANK ( ICNT) 

2  F0RMAT<1X,S13,SX,F8.6,3X,F8.6,7X,F8.6,9X,I2) 

IF( SCORE 2 ( I COUNT) .CT . 0 . *. AND . SCORE 1 ( XCOUNT) . LT. 0 . 7 ) 
♦WRITE  FREE  (10)  *  POSSIRLE  OFF-CENTERED  PARTIAL  TARGET" 
A6  CONTINUE 

WRITE  TO  PRINTER,  VIA  "LPT" 

WRITE  (12,1  A ) MAI NFL ( 1 ) 

WRITE  (12,11) 

WRITE  (12.12) 

1A  FORMATOX, ///////,  30X,  "EVALUATION  RESULTS"  ./.  2SX  ,  "SCENE 
♦FILENAME:  ",S16) 

11  F0RMAT(8X,"  FILE  NAME  8CORE1  SCORE2  TOTAL 

♦SCORE  RANK") 

12  FORMAT  (RX , " -  -  - -  - 

♦ -  - ") 

DO  A7  ICOULT-1 , NUMFILES 
I CNT* I COUNT  ; to  lava  room 

ITEN*ICOUNT* 10 

WRITE (12,3) INFLNM( ITEN) , SCORE! ( ICNT) ,SC0RE2( ICNT) , 

♦SCORE (ICNT), IRANK(ICNT) 

IF( SCORE2 ( ICOUNT) • GT .  0  •  6 . AND • SC0RE1 ( I COUNT) • LT*  0 • 7 ) 
♦WRITE  (12,13) 

13  .FORMATOM  i*  POSSIBLE  OFF-CENTERED  PARTIAL  TARGET") 

3  FORMAT (9X ,S13,3X,F8.t,AX,P8.A,3X,F8.6, 8X , 12 ) 

A 7  CONTINUE 

31  CALL  RESET 

STOP"<7X7X7X7>EVAL  PRINT  LPT  WILL  GIVE  HARD  COPT" 
END 


APPENDI/  D 


Support  Software 

This  Appendix  contains  the  support  software  used 
thoughout  this  thesis.  Included  are  the  following  programs 
and  subroutines: 

VID.FR 

WINDOW. PR 

BI.FR 

HISTO.FR 

NEG.FR 

UNPACK.FR 

repack.fr 

a 


.V  % 
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c********************** **************************************** 

c 

C  VID.FR  -  DC  FORTRAN  4  -  BY  LT.  RICHARD  L.  MILLS  10  AUG  1.964 
C 

C  VID  will  use  the  OCTEK  to  display  files  smaller  then 
C  256  x  256. 

C 

C  RELOAD  COMMAND  (from  DP0:NMXLLS  on  the  NOVA) 

C  RLDR  MI LLS : VI D  MI LLS : VI NDOW  IACF4.LB  NF4PICBUF.LB  FORT. LB 

C 

C  EXECUTION  COMMAND: 

C  VID  (from  NOVA  only) 

C 

c**«* ********** ***************************************** ******* 
C"  PICBUF  STUFF 

PARAMETER  NBUFSZ-300 
PARAMETER  NCOLSMAX-312 
INTEGER  IBUFF3(NBUFSZ) 

INTEGER  I ARRAY (NCOLSMAX ) , I NAME ( 20 ) 

INTEGER  I COL , I ROW , IXP.IYP 
INTEGER  DISP(200),ICT(120) 

INITIALIZE  OCKTEK  DISPLAY  by  Lt.J. Tilley 

IFIELD-0 
IX-0 
IXL-320 
IY-0 
IYL-240 
I ENWV  •«  1  5 
IENVC-15 
ICMAX-11 
CANGLS-0 
GRID-0 

CALL  SXNTRO  ( I CT , 63R , I ER ) 

IF  (IER.EQ. 1 )  GOTO  31 

TYPE  "INTRO  ERROR  CODE:  ",IER 
STOP  "UNABLE  TO  INITIALIZE" 


i 

i  ) 


I 

i 


1  CALL  OPEN( 3 , "IACMON. XB",2,IER) 

IF  (IER.NE.1)  TYPE  "WARNING:  UNABLE  TO  ACCESS  IACMON. XB" 
IF  (IER.EQ. 1)  CALL  LXB  (ICT,3) 

CALL  MPRUN  (ICT.0,1) 

CALL  BOXCUR( ICT^ I XL. IYL) 

CALL  MCTAB ( I CT , I X , I Y ) 

CALL  GREYSCALE  (1CT,1) 

CALL  PXFILL(ICT, 13, 1,320, 1,240) 

End  of  Tllloy'e  inits 


11S 


•V  > 


<fi  u  u  u  o  o  o 


c 

C  ACCEPT  INPUT 
C 

100  ACCEPT"NAME  OF  THE  INPUT  FILE:-->" 

READ ( 1 1 ,40) INAME( 1  ) 

40  FORMAT ( S40 ) 

5  ACCEPT"SIZE ISIZE 

C 

C  TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 
C  (due  to  comtrenti  on  DIRECT  4  INVERSE) 

C 

I F ( I  IZE.EQ. 128)  GO  TO  6 
IF(ISIZE.KQ.64)  CO  TO  6 
IF(ISIZE.EQ.32)  CO  TO  6 
IF( ISIZE. EQ .16)  CO  TO  6 
IF(ISIZE.EQ.B)  CO  TO  6 
CO  TO  5 

CONTINUE  ;  (Ire  Is  okay 
GET  POSITION 
TYPE"  " 

ACCEPT" POS I TI ON  TO  BE  DISPLAYED  (ROW, COLUMN) — >",IYP,IXP 
IXP-IXP+32 

I COL* I S I ZE 
I  ROW" I S I ZE 

PICBUFF  FOR  INPUT  FILE 


NROWS" I S I ZE 
NCOLS-ISIZE 
MBITS-4 
I MODE" 1 

CALL  PFMT( I BUFF 3 , NROWS , NCOLS , MBITS , I MODE) 

CALL,  MAKB(IBUFF3) 

CALL  PICIN(IBUFF3, 1 NAME , 2 HDR ) 

CALL  WINDOV( IBUFF3, ICOL, IROW, IXP , IYP, ICT) 

CALL  RELB( IBUFF3) 

ITOCCLE-O 
TYPE"  " 

ACCEPT" ENTER  1  TO  INCLUDE  ANOTHER  PICTURE:  " , ITOGGLE 
I F( ITOCCLE . EQ . 1 )GO  TO  100 

CALL  RESET 
STOP"<7X7X7>VID" 

END 


vv 


c 

C  VINOOW.FR  -DC  FORTRAN  4-  LT.,J.  TILLEY 
C  UPDATED  BY  LT .  RICHARD  L.  MILLS  10  AUG  19S4 

C 

C  CALED  BY: 

C  VID.FR 

C 

G**«**« ********* ft************************************** ******* 

C  The  subroutine  WINDOW. FR  will  do  the  actual  writing  of 

C  the  video  blocks  to  the  Octek  board. 

c******s*ss*sessee******* *********** ******e**e********eesee*ss 

SUBROUTINE  VINDOVf I BUFFS , ICOL, IROW, IXP, I YP, ICT) 

INTEGER  I  BUFF 3 ( 300),IC0L, I SOW , IXP , IYF  ;cosaon  arrays 

INTEGER  DISP(200) , ICT (120) 

DO  200  J- 1 , IROW 

CALL  GROV( IBUFF3,J,1 ,1C0L,DISP) 

IXPOS-IXP 
I YFOS* I YF+J 

CALL  VVBLK ( ICT , DISP , IXPOS , ICOL , I YPOS , 1 ) 

200  CONTINUE 


o 


RETURN 

END 


o  r»  a 


c 

C  B  I  .  FR  -DC  FORTRAN  5  -  LT  RICHARD  MILLS,  12  S£FT  1984 

C 

C  Th  i  a  prograa  will  rewrite  a  video  file  into  a  binary 

C  for**t  so  that  it  nay  be  displayed  with  the  prograa 

C  PL0T3D . 

C 

C  RELOAD  LINE: 

C  RLDR  BI  UNPACK  0FLIB0 

C 

C  COMMAND  LINE: 

C  BI 

C 

c* ************* ****** *********** ******•***************•***•**< 

INTEGER  TEMP(2S6),IN(1024) 

INTECER  INFLNM( 7 ) , OFLNM2 ( 7 ) 

C 

C  ACCEPT  INPUTS 
C 

ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ(1 1 , l)INFLNM(l) 

1  FOKMAT( S 1 3 ) 

ACCEPT" ENTER  OUTPUT  FILENAME  ->" 

READ( 11,1 )OFLNM2( 1 ) 

CALL  OPEN( 1 , INFLNM, 1 , IER ) 

CALL  CHECK(IER) 

CALL  DFILW  (OFLNM2 , IER ) 

CALL  CF1LW  ( OFLNM2 ,2 , KER ) 

CALL  CNECK(KER) 

OPEN  2,  OFLNM2 , ATT"”OR" , ERR* 1 00 
100  CONTINUE 
5  TYPE"  " 

TYPE"ENTER  SIZE  OF  FILES:" 

ACCEPT”  234,  12B,  44,  OR  32  — >",ISIZE 

TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 

IF(ISIZE.EQ.236)  CO  TO  6 
IPOSIZS.EQ.  128)  CO  YO  4 
IF( ISIZE.ZQ.44)  CO  TO  4 
IP(I3IZE. EQ. 32 }  GO  TO  4 
CO  TO  3 

4  CONTINUE  ;  *is«  It  okay 

I FACTOR" 2 3 4 / 1 S I ZE  ;find  scaling  factor 

ITiMES-64/ I  FACTOR* "2-1 
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TRANSFER  FILE  INTO  NEW  BINARY  FILE 

DO  33  I-O.ITIMES 

CALL  RDBLK< 1 , I, TEMP. I ,IER) 
CALL  CHECK < I ER ) 

CALL  UNPACK(256,TEMP,IN) 

DO  2  J-1,102* 

WRITE  BINARY  (2)  IN(J) 

2  CONTINUE 


33  CONTINUE 


CALL  RESET 
STOP"<7><7><7><7>BI1" 

END 


» 


V  V 


sS 


r-* 
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o  o  o  v*  non 


C  HiSTO.FR  -  DC  FORTRAN  5  -  LT  RICHARD  MILLS,  II  SEPT  1984  * 

C  * 

C  Thia  program  will  exaaine  the  pixel  value*  of  a  video  * 

C  file  by  producing  a  hi itograa  of  the  value*.  If  a  * 

C  Tektronix  4010  graphic*  terminal  i*  available  a  graph  M 

C  of  the  hietogra*  aay  be  constructed.  * 

C  * 

C  RELOAD  LINE:  * 

C  RLDR  HISTO  UNPACK  8CRAPH0  9FLIB?  * 

C  • 

C  COMMAND  LINE:  * 

C  HISTO  * 

C  * 

c* ************************* ********* **************************** a 

INTEGER  TEMP(256) , IN( 1024) 

INTECER  I NFLNM ( 7 ) 

REAL  HI S ( 16) 

DO  40  1-1,16  ; SET  HIS  ARRAY  TO  ZEROS 

HIS(I)-0.0 
40  CONTINUE 

ACCEPT  INPUTS 

ACCEPT "ENTER  FILENAME  ->" 

READ( 1 1 , 1 ) INFLNMf 1 ) 

1  FORMAT ( S 1 3 ) 

CALL  OPEN ( 1 , I NFLNM , 1 , I ER ) 

CALL  CHECK(IER) 

TYPE"  " 

TYPE"ENTER  SIZE  OF  FILE:" 

ACCEPT"  2>4.  128,  44,  OR  32  -->",ISIZE 

TEST  FOR  ONLY  THE  ALLOWABLE  SIZES 

IF(ISIZE.EQ.2S4)  GO  TO  6 
IFCISIZE.EQ. 128)  GO  TO  4 
I F ( ISIZE.EQ.44)  GO  TO  4 
IF< ISIZI.EO.32)  CO  TO  4 
CO  TO  3 

4  CONTINUE  t  *1<*  U  okay 

I  FACTOR- 2 5 4 / I S 1 ZE  (calculate  acalln^  factor* 

ITIMES-44/I FACTOR**  2 - 1 


non 


c 

C  READ  FILE  VALUES 
C 

DO  3  1»0,ITIMES 

CALL  RDSLK< I , I .TEMP, 1 , IER> 

CALL  CHECK ( I ER ) 

CALL  UN PACK ( 256, TEMP, IN) 

DO  2  J-I.102A 

HIS<IN(J)+l)-HIS(IN(J)+l)+I.O 

2  CONTINUE 

3  CONTINUE 
C 

C  TYPE  HISTOGRAM  INFORMATION 
C 

DO  SO  1*1,16 
IOFF-I-1 

TYPE" VALUE : " . IOFF , "  NUMSER : ” ,HI S ( 1  ) 

50  CONTINUE 

ACCEPT" ENTER  1  FOR  A  GRAPH  ”, ITOGGLE 
IF( ITOGGLE.ME. 1 )GO  TO  SI 

PLOT  THE  GRAPH 

CALL  CRPH2 ( "HI STOGRAM  OF  tlXEL  VALUES" , 1 , HI S , U , IE, 1 , YMIN , 
♦YMAX.O) 

ACCEPT  DUMMY  ;KEEP  THE  PROMPT  OFF  SCREEN  FOR  THE  COPYING 

51  CALL  RESET 
STOP 


£*••«****** l •••*** ******* *************************************** 

c 

C  ORIGINAL  PROGRAM:  MEG . PR  BY  R.  SIMMONS 
C  REWRITTEN:  BY  R.  MILLS 

C  TL 1 1>  progtia  Mill  create  the  "negative"  zt  a  video  f  r aae . 

C  Used  for  an  alternate  persective  on  the  hard  copies  and 

C  to,  conserve  lincprinter  ink. 

C 

C  The  cjssand  line  foraat  is: 

C 

C  NEC 

C 

C*******************************»************** ***************** 

DIMENSION  I BLKS ( 16384) , NEC (16) 

ACCEPT  "ENTER  POSITIVE  FILENAME:  " 

READ  (11,5)  POS 
5  FORMAT  (S40) 

ACCEPT  "ENTER  NEGATIVE  FILENAME:  ” 

READ  (11,5)  MEC'M 

OPEN  l.POS  j  .osltlve  file  on  ch.  1 
CALL  DFILW  (M.  ,:ER) 

IF( IER.NE. 1 ) GOTO  90 

90  CALL  CF I LW( NEG ,  3,64, EES)  ;not  exist  on  the  disk 

CALL  CHECK ( KER  ) 

OPEN  2. NEC. ATT* "OC", ERR* 100 

C 

100  CALL  RD1LK(1,0,I BLKS ,64 ,MER)  {Read  64  blocks  into  IBLKS 
IF(MER.NE. 1 )TYPE"  RDBLK  error:”, MER 
TYPE  "WORKING" 

DO  1  J-l, 16384 

1  I BLKS ( J ) •NOT (IBLKS(J))  ;fora  the  negative  iaage 

CALL  WRBLK(2.0, IBLKS, 64, NER)  {Write  the  blocks 
IF(NER.NE. 1 )TYPE"  WRBLK  error:", NER 

C 

TYPE  "FINISHED!" 

CALL  RESET  :Cloae  ell  channele  . 


SUBROUTINE  UNPACK ( N , PIXWORD , PIXELS ) 

C 

C  Written  by  Lt .  Slaaon«(Ref  9)  Version  2 

C 

C  This  subroutine  will  unpack  four  4-bit  integers  from  a 

C  16-bit  integer  word.  The  pixels  in  a  video  file  have  to 

C  be  unpacked  if  each  pixel  is  to  be  operated  on  separately. 

C  ' 

INTEGER  PIXWORD(N) ,PIXELS(<  ,N)  {Four  pixels  per  word 
DO  1  1“ 1 ,N  ; *  N '  allows  higher-order 

DO  1  J-1,4  ;arrays  to  be  passed. 

PIXELS((5-J),I)"15. AND . P IXWORD (I) ; Pi ck  off  right  pixel 
1  P IXWORD ( I  )  ■  I SHFT ( PI XWORD ( I ) ,  -  4  )  {Shift  word  4  bite  ri^ht 

RETURN  {to  pick  off  next  pixel. 

END 


SUBROUTINE  REPACK (N, PIXELS, PXWD ) 


Written  by  Lt .  Sinoni  (Ref  9)  Vereion  2 

This  subroutine  will  repeck  four  A-bit  integer  pixel* 
into  one  16-blt  word  for  uee  by  CHOPS.  Parameter  N 
allows  aore  than  one  A-bit  to  1-word  repacking 
operation  in  each  call  to  REPACK. 

INTEGER  PIXELS(A,N) . PXWD ( N ) 

DO  1  J«1,N  ;Loop  N  times 

PXWD ( J )  »0 
DO  1  '  I»1 , A 

PXWD ( J )■! SHPT( PXWD( J ) , A )  {Shift  pixel  left  in  word 

PXWD( J )«PIXELS(I,J )+PXWD( J )  {then  add  next  pixel  on  right 
RETURN 
END 
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Test  Case  II 


This  test  case  deals  with  locating  the  F-I6*s  in  Figure 

2.2  with  the  norttl,  generalized  F-16  teaplate  of  Figure 
3.7b.  The  Wedging  has  already  been  presented  in  the  text  in 
Figure  2.7.  The  first  pass  correlation,  uaing  the  39x39 
pixel  filter,  was  also  presented  in  Figure  3.11.  Fro*  this 
correlation  output  the  saae  five  potential  targets  of  Figure 

4.2  were  extracted  froa  the  edge  detected  laage,  energy 
noraellzed  and  reccorelated  with  the  teaplate  of  Figure 
3.7b. 

Figure  E.l  shows  the  output  of  the  recorrelation  of 
these  saaller,  energy  noraalised  laages.  The  evaluation  of 
the  peaks  within  these  iaages  are  Included  in  Table  II. 
Note  that  rating  for  the  4  actual  targets  is  rather  low,  but 
it  is  higher  than  the  rating  for  the  iaage  of  the  false 
target. 
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Figure  B.le-c 
a) 


a-e  Second  Fees  Correlation  Output  of  Ceee  II 

a)  Correlation  Output  of  Figure  4.2a  with 

General  F-16  Template  of  Figure  3.7b, 

b)  Correlation  Output  of  Figure  4.2b  with 

General  F-16  Teaplate  of  Figure  3.7b, 

c)  Correlation  Output  of  Figure  4.2c  with 

General  F-16  Teaplate  of  Figure  3.7b., 

(No  reduction  of  FIX  output) 
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Figure  I.ld>«  Second  Feee  Correlation  Output  of  Case  XI. 

d)  Correlation  Output  of  Figure  4. 2d  with 

General  F-16  Teaplate  of  Figure  1.7b, 

e)  Correlation  Output  of  Figure  4.2e  with 

General  F-16  Teaplate  of  Figure  2.7b. 

(No  reduction  of  FIX  output) 
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Table  II 

Evaluation  of  Teat  Case  II 


13 

L- 

V 

| 

g 


FILE  NAME 

FIOEiA.  VD 
FI0E1B.  VD 
FIGE1C.  VD 
FI0E1D.  VD 
FI0E1E.  VD 


EVALUATION  RESULTS 
SCENE  FILENAME:  FIGURE  E.  I 

SCORE 1  8C0RE2  TOTAL  SCORE 


7*6275 

811202 

796275 

827416 

733916 


.  742080 
.783621 
.715959 
. 723781 
.  366479 


. 390900 
.633676 
.  370100 
.  600323 
.  413748 


RANK 

6 

7 

6 

7 

5 


i'N 

:•* 


t c 

t 


K 

% 

N 
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Test  Case  III 


The  leete  of  three  T-34'*,  shown  in  Figure  E.2,  is 
considered  in  this  third  test  cese.  The  edge  extrection  is 
shown  in  Figure  E.3.  Figure  E.4  contains  the  templates  used 
for  the  first  and  the  second  pass  correlation*  The  second, 
saaller  teaplate  was  used  to  keep  the  iaage  sisea  within 
64x64  pixels. 

The  output  of  the  first  pass  correlation  is  given  in 
Figure  E.3  and  the  saaller,  potential  targets  corresponding 
to  the  peaks  are  shown  in  Figure  E.6.  The  second  pass 
correlation  output  is  shown  in  Figure  E.7  and  the  evaluation 
of  these  peaks  is  presented  in  Table  XII. 
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Table  III 


FILE  NAME 


Evaluation  of  Teat  Caae  III 


EVALUATION  RESULTS 
SCENE  FILENAME:  FIGURE  E.  7 

SCORE 1  SCORES  TOTAL  SCORE 


FI9E7A.  VD  .  977903  .  8395 IS  .820964 

FI0E7B.  VD  .  240619  .  731477  .  176007 

POSSIBLE  OFF-CENTERED  PARTIAL  TARGET 
FI0E7C.  VD  .977903  .  710138  .  694446 

FI0E7D.  VD  .  977903  .  494797  .  483364 


t*  *4  W  <0 


T«  a  t  Case  IV 


This  fourth  test  case  exislncs  the  iaage  of  five  F-86's 
shown  in  Figure  E.8.  The  iaage  is  relatively  free  froa 
clutter  except  for  the  contrails  which  show  up  very  well 
after  edge  extraction  as  shown  in  Figure  E.9.  Figure  E.10 
is  the  target  teaplate  used  to  produce  the  first  pass  corre¬ 
lation  output  of  Figure  E.ll.  The  eaall  potential  targets 
are  shown  in  Figure  E.12  and  their  second  pees,  energy 
noraelixed  correlations  are  presented  in  Figure  E.13.  The 
evaluation  of  these  eight  recorrelations  is  given  in  Table 
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Figure  E.  10  F-»*  Template. 
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Figure  I. 11  Firat  Pati  Correlation  Output  of  Taat  Case  IV. 


Fl|ur«  t.l2e-h  Snell  Piles  Corresponding  to  Peeks  in 

Figure  K.ll.  (No  reduction  of  FIX  output) 


Table  IV 


Evaluation  of  Taat  Caaa  IV 


FILE  NAME 

FI0E13A.  VO 
FI0E13B.  VD 
FICE13C.  VD 
F10E13D.  VD 
POSSIBLE 
FI0E13E.  VD 
FIGE13F.  VD 
FIDE 130.  VD 
FI0E13H.  VD 


EVALUATION  RESULTS 
SCENE  FILENAME:  FIGURE  E.  13 

SCORE 1  8C0RE2  TOTAL  SCORE 


RANK 


.  *77903 
.  *77903 
.  000000 
.  oooooo 


.  89*001 
.  670238 
.  961368 
.  696622 


OFF-CENTERED  PARTIAL  TARGET 
.977903  .  963494 

.  449G06  .  990302 

.  139348  .  490834 

.  977903  .  489838 


879136 

699427 

OOOOOO 

oooooo 

991003 

249329 

068396 

479014 


_ 
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Test  Css*  V 


This  fifth  test  esse  exseines  the  image  of  four  1-52's 


shown  in  Figure  E.14.  The  iaagi 


itfically  crested 


using  s  scene  generating  progrse  by  Lt  J.  Tilley,  which 
superieposed  various  aircraft  images  onto  a  noisy  cloud 
background.  The  edge  extraction  of  the  scene  is  shown  in 


Figure  E. li.  Figure  E.16  is  the  edged  target  template  used 
to  produce  the  first  pass  correlation  output  of  Figure  E.17. 


The  small  potential  targets,  corresponding  to  the  correla¬ 
tion  peaks,  are  shown  in  Figure  E.1S.  The  second  pass, 
energy  normalised  correlations  are  shown  in  Figure  E.19. 
The  evaluation  of  these  four  recorrelations  is  presented  in 
Table  V. 
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Flyura  1.18  Small  Piles  Corrssi 

Figure  E.17.  I 
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Figure  E.19  Second  Pees  Correlation  Output  of  Teat  Caee  V. 

(No  reduction  of  FIX  output) 
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Table  V 

Evaluation  of  Teat  Case  V 


FILE  NAME 


EVALUATION  RESULTS 
SCENE  FILENAME:  FIGURE  E.  19 

SCORE 1  SC0RE2  TOTAL  SCORE 


RANK 


FI0E19A.  VD 
FI0E19B.  VD 
FIGE19C.  VD 
FIGE19D.  VD 


889519 

889515 

889515 

889515 


.  978030 
.  920735 
.984499 
.  970009 


.  869972 
.819008 
.875726 
.862837 
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Teat  Cass  VI 


The  lsags  of  two  B-32'e.  and  two  KC-135's  on  a  cluttered 
backgr  ound,  shown  in  Figure  E.20,  is  exasined  in  Test  Case 
VI.  This  xaage  was  also  arltflcally  produced.  The  edge 
extraction  of  the  scene  is  shown  in  Figure  E.21.  Note  that 
the  I-J2  targets  were  hidden  in  the  clouds  while  the  KC-135 
decoys  were  placed  in  relatively  clear  view.  Although  this 
thesis  doss  not  eaphasise  the  ability  to  distinguish 
the  difference  between  different  types  of  targets,  this  test 
case  was  used  to  see  what  the  algoritha's  perforaance  Was  in 
such  a  situation. 

The  edged  target  teaplate  used  in, this  test  case  was 
the  saae  B-52  teaplate  in  Figure  E.I6.  The  first  pass 
correlation  output  is  shown  in  Figure  E.22.  The  saall 
potential  targets,  corresponding  to  the  correlation  peaks, 
are  shown  in  Figure  E.23.  The  second  pass,  energy  noraa- 
llzed  correlations  are  shown  in  Figure  E.24.  The  evaluation 
of  these  four  rdcorrelav *«%••  as  prsssntsd  in  Tsbls  VI.  Note 
thet  the  laages  of  the  hidden  ••32* s  did  rats  higher  the  the 
lasges  of  the  KC-13S'e  that  were  in  clear  view. 
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Table  VI 


Evaluation  of  Taat  Case  VI 


FILE  NAME 


EVALUATION  RESULTS 
SCENE  FILENAME:  FIGURE  E.  24 

SCORE 1  SC0RE2  TOTAL  SCORE 


RANK 


FI0E24A.  VD 
FI0E24B.  VD 
FI0E24C.  VD 
FI0E24D.  VD 


977903 

977903 

977903 

977903 


.  873436 
.  655509 
. 980772 
. 720717 


.  854136 
.641020 
.959100 
.704791 


Test  Case  VII 


The  final  test  case  Included  In  this  thesis  la  of  an 
inage  that  does  not  contain  any  targets.  This  is  the 
cluttered  cloud  background,  used  in  the  two  previous  test 
cases,  shown  in  Figure  E.2S.  The  edge  extraction  of  the 
scene  is  shown  in  Figure  E.26. 

The  edged  target  tenplate  used  in  this  test  case  was 
the  sase  1-52  teaplate  in  Figure  E.16.  The  first  pass 
correlation  output  is  shown  in  Figure  E.27  The  saall  poten¬ 
tial  targets,  corresponding  to  the  weak  correlation  peaks, 
are  shown  in  Figure  E.28.  Four  peaks  were  chosen  as  being 
locations  of  possible  targets.  The  second  pass,  energy 
noraallted  correlations  are  shown  in  Figure  E.29.  The 
evaluation  of  these  four  recorrelations  is  given  in  Table 
VI.  Note  that  two  areas  of  noise  ranked  higher  than  the 
other  two  areas  of  noise  due  to  their  isolated  bnergy 
distribution  in  the  center  of  the  scene  which  notched  with 
the  central  energy  of  the  aircraft  tenplate. 
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Figure  E.2S 


*••11  Piles  Corresponding  to  Peaks  in 

Figure  K.27.  (No  reduction  of  FIX  output) 
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Figure  E. 29  Second  Pan  Correlation  Output  of  Teat  Cate  VII, 

, (No  reduction  of  FIX  output) 


FILE  NAME 

FICE29A.  VD 
FI0E29B.  VD 
FIGE29C.  VD 
FI0E29D.  VD 


Table  VII 

Evaluation  of  Teat  Caaa  VII 


EVALUATION  RESULTS 


SCENE  FILENAME:  FIGURE 
SCORE 1  SC0RE2 

.125838  .416693 
.  950589  .  560848 
.  S58510  .649143 
.  276173  .  457362 


E.  29 

TOTAL  SCORE  RANK 

.  052436  1 

.  533136  6 

.  557295  6 

.  126311  2 
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